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Abstract 

A set of boundary conditions defining an undistorted, non-rotating isolated horizon 
are specified in general relativity. A space-time representing a black hole which is itself 
in equilibrium but whose exterior contains radiation admits such a horizon. However, 
■ the definition is applicable in a more general context, such as cosmological horizons. 

Q\ . Physically motivated, (quasi-)local definitions of the mass and surface gravity of an 

*q \ isolated horizon are introduced and their properties analyzed. Although their defini- 

tions do not refer to infinity, these quantities assume their standard values in the static 
Sh . black hole solutions. Finally, using these definitions, the zeroth and first laws of black 

hole mechanics are established for isolated horizons. 

> ■ 

• rH 

h : 1 Introduction 

a i 

The similarity between the laws of black hole mechanics and those of ordinary thermodynam- 
ics is one of the most remarkable results to emerge from classical general relativity |], @, |3|. Q • 
However, the original framework was somewhat restricted and it is of considerable interest 
to extend it in certain directions, motivated by physical considerations. The purpose of this 
paper is to present one such extension. 

The zeroth and first laws refer to equilibrium situations and small departures therefrom. 
Therefore, in this context, it is natural to focus on isolated black holes. In the standard treat- 
ments, these are generally represented by stationary solutions of field equations, i.e, solutions 
which admit a time-translational Killing vector field everywhere, not just in a small neigh- 
borhood of the black hole. While this simple idealization is a natural starting point, it seems 
to be overly restrictive. Physically, it should be sufficient to impose boundary conditions at 
the horizon which ensure only that the black hole itself is isolated. That is, it should suffice 
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Figure 1: (a) A typical gravitational collapse. The portion A of the horizon at late times 
is isolated. The space-time M. of interest is the triangular region bounded by A, and 
a partial Cauchy slice M. (b) Space-time diagram of a black hole which is initially in 
equilibrium, absorbs a small amount of radiation, and again settles down to equilibrium. 
Portions Ai and A2 of the horizon are isolated. 



to demand only that the intrinsic geometry of the horizon be time independent, whereas the 
geometry outside may be dynamical and admit gravitational and other radiation. Indeed, 
we adopt a similar viewpoint in ordinary thermodynamics; in the standard description of 
equilibrium configurations of systems such as a classical gas, one usually assumes that only 
the system under consideration is in equilibrium and stationary, not the whole world. For 
black holes in realistic situations, one is typically interested in the final stages of collapse 
where the black hole is formed and has 'settled down' or in situations in which an already 
formed black hole is isolated for the duration of the experiment (see figure 1). In such situ- 
ations, there is likely to be gravitational radiation and non-stationary matter far away from 
the black hole, whence the space-time as a whole is not expected to be stationary. Surely, 
black hole mechanics should incorporate such situations. 

A second limitation of the standard framework lies in its dependence on event horizons 
which can only be constructed retroactively, after knowing the complete evolution of space- 
time. Consider for example, figure 2 in which a spherical star of mass M undergoes a 
gravitational collapse. The singularity is hidden inside the null surface Ai at r = 2M which 
is foliated by a family of marginally trapped surfaces and would be a part of the event horizon 
if nothing further happens. Suppose instead, after a very long time, a thin spherical shell of 
mass 5M collapses. Then Ai would not be a part of the event horizon which would actually 
lie slightly outside Ai and coincide with the surface r = 2(M + 5M) in the distant future. 
On physical grounds, it seems unreasonable to exclude Ai a priori from thermodynamical 
considerations. Surely one should be able to establish the standard laws of mechanics not 
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Figure 2: A spherical star of mass M undergoes collapse. Later, a spherical shell of mass 
SM falls into the resulting black hole. While Aj and A 2 are both isolated horizons, only A 2 
is part of the event horizon. 



only for the event horizon but also for A x . 

Another example is provided by cosmological horizons in de Sitter space-time ||. In this 
case, there are no singularities or black-hole event horizons. On the other hand, semi-classical 
considerations enable one to assign entropy and temperature to these horizons as well. This 
suggests the notion of event horizons is too restrictive for thermodynamical analogies. We 
will see that this is indeed the far as equilibrium properties are concerned, the 

notion of event horizons can be replaced by a more general, quasi-local notion of 'isolated 
horizons' for which the familiar laws continue to hold. The surface Ai in figure 2 as well as 
the cosmological horizons in de Sitter space-times are examples of isolated horizons. 

In addition to overcoming these two limitations, the framework presented here provides 
a natural point of departure for quantization and entropy calculations || [7|, ||. In con- 
trast, standard treatments of black hole mechanics are often based on differential geometric 
identities and are not well-suited to quantization. 

At first sight, it may appear that only a small extension of the standard framework 
based on stationary event horizons is needed to overcome the limitations discussed above. 
However, this is not the case. For example, in the stationary context, one identifies the 
black-hole mass with the ADM mass defined at spatial infinity. In the presence of radiation, 
this simple strategy is no longer viable since even radiation fields well outside the horizon 
also contribute to the ADM mass. Hence, to formulate the first law, a new definition of 
the black hole mass is needed. Similarly, in the absence of a globally defined Killing field, 
we need to generalize the notion of surface gravity in a non-trivial fashion. Indeed, even 
if space-time happens to be static in a neighborhood of the horizon — already a stronger 
condition than contemplated above — the notion of surface gravity is ambiguous because the 
standard expression fails to be invariant under constant rescalings of the Killing field. When 
a global Killing field exists, the ambiguity is removed by requiring the Killing field be unit at 
infinity. Thus, contrary to intuitive expectation, the standard notion of the surface gravity of 
a stationary black hole refers not just to the structure at the horizon, but also to infinity. This 
'normalization problem' in the definition of the surface gravity seems especially difficult in the 
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case of cosmological horizons in (Lorentzian) space-times whose Cauchy surfaces are compact. 
Apart from these conceptual problems, a host of technical issues must also be resolved. In 
Einstein-Maxwell theory, the space of stationary black hole solutions is finite-dimensional 
whereas the space of solutions admitting isolated horizons is m/imte-dimensional since these 
solutions also admit radiation. As a result, the introduction of a well-defined action principle 
is subtle and the Hamiltonian framework acquires qualitatively new features. 

The organization of this paper is as follows. In Section ^| we recall the formulation of 
general relativity in terms of SL(2, C)-spin soldering forms and self-dual connections for 
asymptotically flat space-times without internal boundaries. In Section |3], we specify the 
boundary conditions which define non-rotating isolated horizons and discuss several exam- 
ples. The primary focus in Section 3, as in the rest of the paper, is on Einstein-Maxwell 
theory, although more general matter is also considered. The consequences of these bound- 
ary conditions which are needed to obtain the laws governing isolated horizons are discussed 
in Section |j. Using this structure, we introduce in Section ^ the notion of the surface gravity 
k of an isolated horizon without any reference to infinity and prove the zeroth law. 

The action principle and the Hamiltonian formulation are discussed in Section |6|. The 
Hamiltonian has a bulk term and two surface terms, one at infinity and the other at the 
isolated horizon A. The bulk term is the standard linear combination of constraints and the 
surface term at infinity yields the ADM energy. In Section [7], we argue the horizon surface 
term in the Hamiltonian should be identified with the mass Ma of the isolated horizon. In 
particular, in the situation depicted in figure 1(a) we show Ma is the difference between 
the ADM energy and the energy radiated away through all of null infinity (provided certain 
assumptions on the structure at i + hold). That is, although the definition of Ma uses only 
the structure available at the horizon, it equals the future limit of the Bondi mass defined 
entirely on ^ + in the case under consideration. This is precisely what one might expect on 
physical grounds in the presence of radiation. Finally, having k and Ma at our disposal, we 
establish the first law in Section 

The overall viewpoint and the boundary conditions of Section |3] are closely related to 
those introduced by Hayward in an interesting series of papers [[J which also aims at provid- 
ing a more physical framework for discussing black holes. There is also an inevitable overlap 
between this paper and @ in which a Hamiltonian framework (in terms of real SU(2) con- 
nections) is constructed with an eye towards quantization. However, by and large the results 
presented here are complementary to those obtained in 0. Even when there is an overlap, 
the material is presented from a different angle. The relation between these three papers is 
discussed in Section |^. Technical, background material needed at various stages is collected 
in the three appendices. 



A brief summary of our main results was given in |]T0 
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2 Mathematical Preliminaries 



In this paper we will use the formulation of general relativity as a dynamical theory of 
connections |TT| rather than metrics. Classically, the theory is unchanged by the shift to 



connection variables. However, at present the connection variables appear to be indispens- 
able for non-perturbative quantization [|12|] . In particular, the entropy calculation for isolated 
horizons has been carried out only within this framework ||. Therefore, for uniformity, we 
will use the connection variables in the main body of this paper. However, as indicated in 
Section [|, all these results can also be derived in the framework of tetrad dynamics. 

In order to fix notation and conventions and to acquaint the reader with the basics of con- 
nection dynamics, we will provide here a brief review of this formulation in an asymptotically 
flat space-time without interior boundary. (The modifications required to accommodate a 
non-zero cosmological constant are discussed at the end.) For further details see, e.g., |I2 . 



2.1 Connection Variables 

Fix a four- dimensional manifold M.. In Einstein-Maxwell theory, the basic fields consist of 
the triplet of asymptotically flat, smooth fields (a a AA ' , A a AB , A a ). Here, o a AA ' is an anti- 
Hermitian soldering form for SL(2, C) spinors, A a AB is a self dual SL(2, <C) connection acting 
only on unprimed spinors and A a is the U(l) electro-magnetic connection. The action for 
the theory is given by 

S(M,A) = ^ / Tr[£AF] + i- / F A *F + -J- / Tr[E A A]. (2.1) 

07TG J M 07T J M 07rG JCoo 

The 2-form fields £ are given by Y, AB = a AA ' A a A , B , F is the curvature of the connection A, 
F is the electro-magnetic field strength, and is the time-like cylinder at spatial infinity. 
We define a metric g ab of signature — h ++ on this manifold via g a b = u a AA c aa'- With 
respect to this metric, the 2-form fields £ are self dual. (For details, see |13|.) 

Let us consider the equations of motion arising from this action. Varying the action with 
respect to the connection A, one obtains 

£>£ = 0. (2.2) 

This implies the connection D defined by A has the same action on unprimed spinors as the 
self dual part of the connection V compatible with the soldering form, V a & AA = 0- When 
this equation is satisfied, the curvature F is related to the Riemann curvature by: 

TP AB — v AB f? cd (1 i\ 

fab - ~~^cd K ab ■ \ Z - 6 ) 

Varying the action Q2.1D with respect to a and taking into account the compatibility of A 
with a, we obtain a second equation of motion 

G ab = 8nGT ab (2.4) 
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where G a b is the Einstein tensor of g a b and T a b is the standard stress-energy tensor of the 
Maxwell field F@,§. 

Next, let us consider the equations of motion for the Maxwell field, 

d¥ = 0, and cfF = 0. (2.5) 

Since F is the curvature of the U(l) connection A, the first Maxwell equation d¥ = is 
an identity. If one varies equation (|2.1| ) with respect to A, one obtains the second Maxwell 
equation, d*F = 0. Thus, the equations of motion which follow from the action ( |2.1[ ) are the 
same as those given by the usual Einstein-Hilbert-Maxwell action; the two classical theories 
are equivalent. 



2.2 Hamiltonian Formulation 

To pass to the Hamiltonian description of the theory, it is necessary to re-express the action 
in terms of 3-dimensional fields. Let us assume the space-time M. is topologically M x R. 
Introduce a 'time function' t which agrees with a standard time coordinate defined by the 
asymptotically Minkowskian metric at infinity. A typical constant t leaf of the foliation will 
be denoted M. Fix a smooth time-like vector field t a , transverse to the leaves M such that: 
(i) £ a V a t = 1 and (ii) t a tends to the unit time translation at spatial infinity. We will denote 
the future directed, unit normal to the leaves M by r a . The intrinsic metric on the 3-surfaces 
M is g a b := 4 g ab + T a r fc .[] As usual, by projecting t a into and orthogonal to M, we obtain the 
the lapse and shift fields, N and N a respectively: t a = Nr a + N a . 

We are now in a position to define the basic phase space variables of the theory. They 
are simply the pull backs to the space-like 3-surfaces M of the space-time variables A, £ and 
A, together with the electric field two-form E. To perform the Legendre transform, note the 
pull back of the soldering form a to M induces an SU(2) soldering form and the pull back 
of the connection A induces a complex-valued, SU(2) connection on spatial spinors. More 
precisely, we have 

o- AB o,/oj> 4 o- AA ' rr B > ^ 4_ AA! • ft „ A BA' ^ AA' 

4 A A B = A 0A B dt + A A B and 4 A = A dt + A, 

where t aa ' := 4 a a AA 'r a , r a4 A a AB = and r a 4 A a = 0. (Note that this decomposition of 
connection 1-forms uses the space-time metric.) Using these definitions, one arrives at the 
following 3+1 decompositions of S and the gravitational and electro-magnetic field strengths: 

= X — %N\p2o A dt 
4 F = {-A + D{tM) + Nj F) Adt + F (2.7) 
4 F = (-A + d{t A A) +NjdA)Adt + dA, 



1 In the discussion of the Legendre transform and the Hamiltonian, both in this section and section ^, 
the four dimensional fields will carry a superscript 4 preceding the field and all other fields will be assumed 
to be three-dimensional, living on the space-like surface M. 
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where N is the shift field. The electric field E and magnetic field B are defined as usual to 
be the pull backs to the space-like hypersurface M of *F and F respectively. 

The phase space of the theory consists of quadruples (A a AB , Y> ab AB , A a , E^). These fields 



satisfy the standard falloff conditions [T^|. To specify them, let us fix an SU(2) soldering 
form (jonM such that the 3-metric g a t = P a AB &bAB is flat outside of a compact set. Then 
the quadruple of fields is required to satisfy: 

M(6,(j))\ o 




Ti- 



er 



a A r , 



O 



o 



o 



1 



A=Q — 



E = O — 



(2.8) 



where r is the radial coordinate defined by the flat metric g a b- 

The action can be re-expressed in terms of these fields using equations (|2.6| ) and (|2.7|) . 
From this action, it is straightforward to read off the Hamiltonian and symplectic structure. 



H, 



m 8nG 

+ 
1 

+ 



Tr[(t. 4 A)DS] + 
N 



47T 

(Ti[V2a A F] 



(t 4 A)dE + 



;SA(iVj F) 



8nG 



SnG 

I Tr[V2Na AA + i(N j A)E] 

•J Snr> 



8nG 

G(E A *E + dA A *dA)) 



4tt 



E A (N j dA) 



8ttG 



Tt^A A 5 2 S - 5 2 A A M + — 

At 4:71 



M 



SxA A 5 2 E - 5 2 A A 5xE 



(2.9) 

As always in general relativity, the Hamiltonian takes the form of constraints plus boundary 
terms. The constraints consist of two Gauss law equations, one for the self dual two form E 
and the other for the electric field E, together with the standard vector and scalar constraints. 
When the constraint equations are satisfied, the term at infinity equals t a P a where P a is the 
ADM 4-momentum. In our signature, t a P a is negative, so —t a P a = E ADM , the standard 
ADM energy. The equations of motion are just Hamilton's equations: 



6H = Cl(6, X s ) 



(2.10) 



These are the field equations fl2.4p and fl2.5| ) in a 3+1 form, expressed in terms of the canonical 
variables. 

To conclude, we note two modifications which occur if the cosmological constant is non- 
zero. First, there is an extra term proportional to A Tr[S A E] in the action (|2.1|) , where A 
is the cosmological constant. This contributes a term proportional to ATra A a A a in the 
bulk term of the Hamiltonian ( |2.9| ). Second, the boundary conditions are modified. If A is 
positive, it is natural to assume the Cauchy surfaces are compact, whence there are no falloff 
conditions or boundary terms in any of the expressions. If A is negative, the dynamical fields 



approach asymptotically their values in the anti-de Sitter space-time [15 . 
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3 Boundary Conditions 



In this section, we specify the boundary conditions which define an undistorted, non-rotating 
isolated horizon. As explained in the Introduction, the purpose of these boundary conditions 
is to capture the essential features of a non-rotating, isolated black hole in terms of the 
intrinsic structure available at the horizon, without any reference to infinity or to a static 
Killing field in space-time. However, the boundary conditions model a larger variety of 
situations. For example, the cosmological horizons in de Sitter space-times || are isolated 
horizons, even though there is no sense in which they describe a black hole. As a result, the 
usual mechanics of cosmological horizons will be reproduced within the framework of isolated 
horizons. Other examples involve space-times admitting gravitational and electro-magnetic 



radiation such as those described in Section 3.2 



The physical situation we wish to model is illustrated by the example of figure [TJ(a). 
The late stages of the collapse pictured here should describe a non-dynamical, isolated black 
hole. However, a realistic collapse will generate gravitational radiation which must either be 
scattered back into the black hole or radiated to infinity. Physically, one expects most of 
the back-scattered radiation will be absorbed rather quickly and, in the absence of outside 
perturbations, the black hole will 'settle down' to a steady-state configuration. This picture 
is supported by numerical simulations. The continued presence of radiation elsewhere in 
space-time however implies M. cannot be stationary. As a result, the usual formulations of 
black hole mechanics in terms of stationary solutions cannot be easily applied to this type 
of physical black hole. Nevertheless, since the portion A of the horizon describes an isolated 
black hole one would hope to be able to formulate the laws of black hole mechanics in this 
context. We will see in Sections 5 and 8 that this is indeed the case. 



3.1 Definition 

We are now in a position to state our boundary conditions. Since this paper is concerned 
primarily with the mechanics of isolated horizons in Einstein-Maxwell space-times, condi- 
tions on gravitational and Maxwell fields are specified first and more general matter fields 
are treated afterwards. 

A non-rotating isolated horizon is a sub-manifold A of space-time at which the following 
five conditions hold:Q 

(I) A is tope-logically S 2 xR and comes equipped with a preferred foliation by 2-spheres 
Sa and a ruling by lines transverse to those 2-spheres. 

These preferred structures give rise to a 1-form direction field [n a ] and a (future- 
directed) vector direction field [£ a ] on A. Furthermore, any n a £ [n a ] is normal to a 

2 Throughout this paper, = will denote equality at points of A. For fields defined throughout space-time, 
a single left arrow below an index will indicate the pull-back of that index to A, and a double arrow will 
indicate the pull-back of that index to the preferred 2-sphere cross-sections Sa of A introduced in condition 
H For brevity of presentation, slightly stronger conditions were used in |ic|j . 
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foliation of A by 2-spheres and we further impose dn = 0. As a result, the equivalence 
class [n a ] is defined with respect to rescaling only by functions which are constant on 
each Sa- A function with this property will be said to be spherically symmetric. We 
'tie together' the normalizations of the two direction fields by fixing £ a n a = — 1 with 
£a ^ |-£oj_ ^nis i ea ves a single equivalence class [£ a ,n a ] of direction fields subject to the 
relation 

(e a ,n a )^(F- 1 r,Fn a ), (3.1) 
where F is any positive, spherically symmetric function on A. 

(II) The soldering form <J a AA ' gives rise to a metric in which A is a null surface with [£ a ] 
as its null normal. 

So far the 1-forms n a are defined intrinsically on the 3- manifold A. They can be 
extended uniquely to 4-dimensional, space-time 1-forms (still defined only at points of 
A) by requiring that they be null. We will do so and, for notational simplicity, denote 
the extension also by n a . Given any null pair (£ a ,n a ) satisfying £ a n a = —1, it is easy 
to show there exists a spin basis (i A , o A ) satisfying |£ A oa| — 1 such that 

t = ia a AA , A o A ' and n a = ia AA 'i A l A >. (3.2) 

We will work by fixing a spin dyad (i A , o A ) with l a oa — +1 once and for all and 
regarding ( |3.2| ) as a condition on the field o a AA . Finally, using this spin dyad, we can 
complete (£ a , n") to a null tetrad with the vectors 

m a = i a a AA ,o A l A ' and rh a = ia a AA ,L A d A ' (3.3) 

tangential to the 2-spheres 

(III) The derivatives of the spin dyad are constrained by 

o A VaO A = and l A V = fJ>fh a , (3.4) 

where \i is a real, nowhere vanishing, spherically symmetric function, and V a denotes 
the unique torsion-free connection compatible with u a AA . The function \i is one of the 
standard Newman- Penrose spin coefficients. 

(IV) All equations of motion hold at A. In particular: 

IVa. The SL(2, <D) connection is compatible with the soldering form: D a X A = VqA^ 4 . 

3 Since the topology of A is S 2 x H, the tetrad vectors m a and fh a — and the spin-dyad (i A , o A ) — fail 
to be globally defined. Thus, when we refer to a fixed spin-dyad, we mean dyads which are fixed on two 
patches and related by a gauge transformation on the overlap. Our loose terminology is analogous to the 
one habitually used for (spherical) coordinates on a 2-sphere. 
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IVb. The Einstein equations hold: G a b + ^9ab — 8irGT a b, where T a b is the stress-energy 
tensor of the matter fields under consideration and G a b is the Einstein tensor of 
the metric compatible with o~ a AA . 

IVc. The electro-magnetic field strength F satisfies the Maxwell equations: dF = 
and cf F = 0. 

(V) The Maxwell field strength F has the property that 



1 m a fh b (F-i*F) ab (3.5) 



is spherically symmetric. 

These five conditions define a non-rotating isolated horizon. Each is imposed only at 
the points of A. Let us now discuss the geometrical and physical motivations behind the 
boundary conditions, and see how they capture the intuitive picture outlined above. 

Conditions | and [TJ are straightforward. The first is primarily topological and fixes the 
kinematical structure of the horizon. (While the S 2 xR topology is the most interesting one 
from physical considerations, most of our results go through if S 2 is replaced by a compact, 
2- manifold of higher genus. This issue is briefly discussed at the end of Section [5].) The 
meaning of the preferred foliation of A will be made clear in the discussion of condition [Tl I 
below. The meaning of the preferred ruling can be seen immediately in condition |lj which 
simply requires that A be a null surface, and [£ a ] its null generator. Thus, the preferred 
ruling singles out the null generators of the horizon. 

Condition [iy| is a fairly generic dynamical condition, completely analogous to the one 
usually imposed at null infinity: Any set of boundary conditions must be consistent with the 
equations of the motion at the horizon. It is likely that this condition can be weakened, e.g., 
by requiring only that the pull-backs to A of the equations of motion should hold. However, 
care would be needed to specify the precise form of equations which are to be pulled-back 
since pull-backs of two equivalent sets of equations can be inequivalent. We chose simply to 
avoid this complication. 

Conditions |, |TJ and [IV] are weak; in particular, they are satisfied by a variety of null 



surfaces in any solution to the field equations. It is condition |T| which endows A with 
the structure of an isolated horizon. Technically, this condition restricts the pull-back of 
the self-dual connection compatible with o~ a . (Note the pull-backs to A are important 
because we have introduced the dyad (i A , o A ) only at the surface itself.) Geometrically, it is 
equivalent to requiring the pairs (£ a ,n a ) G [£ a ,n a ] to have the following properties: 

1. l a is geodesic, twist-free, shear-free and divergence- free. 

2. n a is twist-free, shear-free, has nowhere vanishing, spherically symmetric expansion 

Q(n) = 2/i- (3.6) 
and vanishing Newman- Penrose spin coefficient 7r := l a m b 'V a nb on A. 
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As we will see in the next section, the only independent consequence of condition |T| for 
i a is its vanishing divergence; the rest then follows from simple geometry. The vanishing 
divergence of l a is equivalent to the vanishing expansion of the horizon, We will see in Section 



Pjthat this implies there is no flux of matter falling across A, which in turn captures the idea 
that the horizon is isolated. Collectively, the consequences of condition |TJ for n a imply the 
horizon is non-rotating and its intrinsic geometry is undistorted. Requiring \i to be nowhere 
vanishing is equivalent to requiring the expansion of n a to be nowhere vanishing. In black- 
hole space-times, we expect this expansion to be strictly negative on future horizons and 
strictly positive on past horizons. On cosmological horizons, both signs are permissible. In 
section [| we will impose further restrictions tailored to different physical situations. Finally, 
it is easy to verify that these conditions on n a suffice to single out the preferred foliation 
of A. Thus, we could have just required the existence of a foliation satisfying the first three 



conditions, and used [/7| to conclude the foliation is unique. 

Next, let us discuss the condition [V] on the Maxwell field. At first sight, this requirement 
seems to be a severe restriction. However, if the Newman-Penrose component 0o of the 
electro-magnetic field, representing 'the radiation field traversing A', vanishes in a neighbor- 
hood of A and boundary conditions I through IV hold, then <f>\ automatically satisfies |V} 
Heuristically, this feature can be understood as follows. From the definition of <fti in (|3.5|) , 
one can see condition requires the electric and magnetic fluxes across the horizon to be 
spherically symmetric. If this condition did not hold, one would intuitively expect a non- 
rotating black hole to radiate away the asphericities in its electro-magnetic hair, giving rise 
to a non- vanishing 0o- Therefore, it is reasonable to expect |V| would hold once the horizon 
reaches equilibrium. 

Let us now consider the generalization of the boundary conditions to other forms of mat- 
ter. Conditions [lVc| and |V| refer to the Maxwell field; the rest involve only the gravitational 



degrees of freedom and are independent of the matter fields present at the horizon. We will 
now indicate how these two conditions must be modified in the presence of more general 
forms of matter. (For a discussion of dilatonic couplings, see 0, 0.) First, note that con- 
dition |iy] is unambiguous. Hence, [IVcj would simply be replaced with the field equations of 



the relevant matter. Condition [V], on the other hand, is more subtle. In fact, there is no 
universal analog of (|3.5|) which applies to arbitrary matter fields; the boundary conditions 
used in place of condition |V| may vary from case to case. There are, however, two general 
properties which any candidate matter field and its associated boundary conditions must 
possess: 

(V) For any (£ a ,n a ) G [£ a ,n a ], a matter field must satisfy 
V'a. The stress-energy tensor is such that 

r = -T\t (3.7) 
is causal, i.e., future-directed, time-like or null. 
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V'b. The quantity 

e := T ab tn b (3.8) 

is spherically symmetric on A. 
The first of these properties, |V'a|, is an immediate consequence of the (much stronger) 



dominant energy condition which demands — T a b k b be causal for any causal vector k b . Like 
any energy condition, this is a restriction on the types of matter which may be present near 
the horizon. On the other hand, the property |V'b| is a restriction on the possible boundary 
conditions which may be imposed on matter fields at the horizon. 

We conclude this section with a remark on generalizations of these boundary conditions. 
Although our framework is geared to the undistorted, non-rotating case, only the require- 
ments on n a in condition [TTT] and the symmetry condition on the Maxwell field would have 
to be weakened to accommodate distortion and rotation [T/J. Specifically, it appears that, in 
presence of distortion, ri will not be spherically symmetric and in presence of rotation n will 
not vanish. However, it appears that a more general procedure discovered by Lewandowski 
will enable one to introduce a preferred foliation of A and naturally extend the present 
framework to allow for distortion and rotation. 

Finally, in light of the sphericity conditions on 9r n ) and <p±, one may be tempted to call 
our isolated horizons 'spherical'. However, the Newman-Penrose curvature components ^3 
and \l/4 and the Maxwell field component 02 need not be spherical on A for our boundary 
conditions to be satisfied. Therefore, the adjectives 'undistorted and non-rotating' appear 
to be better suited to characterize our isolated horizons. 



3.2 Examples 

It is easy to check that all of these boundary conditions hold at the event horizons of Reissner- 
Nordstrom black hole solutions (with or without a cosmological constant). Similarly, they 
hold at cosmological horizons in de Sitter space-time. Furthermore, if one considers a spher- 
ical collapse, as in figure 2, they hold both on A x and A 2 (at suitably late times). 

In the non- rotating context, these cases already include situations normally considered 
in connection with black hole thermodynamics. However, the isolated horizons in these 
examples are also Killing horizons for globally defined, static Killing fields. We will now 
indicate how one can construct more general isolated horizons. 

First, an infinite-dimensional space of examples can be constructed using Friedrich's 



results pl|, and Rendall's extension [[LP thereof, on the null initial value formulation (see 
figure 3). In this framework, one considers two null hypersurfaces A and M with normals 
i a and n a respectively, which intersect in a 2-sphere S. (At the end of the construction A 



will turn out to be a non-rotating isolated horizon.) In a suitable choice of gauge ||18|| , the 
free data for vacuum Einstein's equations consists of ^0 011 A; on A/"; and, the Newman 
Penrose coefficients A, cr, it, Re [jjl] , Re [p] as well as the intrinsic metric 2 g ab on the 2-sphere 
S. Given these fields, there is a unique solution (modulo diffeomorphisms) to the vacuum 
Einstein's equations in a neighborhood of S bounded by (and including) the appropriate 
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Figure 3: Space-times with isolated horizons can be constructed by solving the characteristic 
initial value problem on two intersecting null surfaces, A and M . The final solution admits 
A as an isolated horizon. Generically, there is radiation arbitrarily close to A and no Killing 
fields in any neighborhood of A. 




Figure 4: A space-time Ai with an isolated horizon A as internal boundary and radiation 
field in the exterior can be obtained by starting with an asymptotic region of Kruskal space- 
time and modifying the initial data on the partial Cauchy surface M. While the new metric 
continues to be isometric with the Schwarzschild metric in a neighborhood of A, it admits 
radiation in a neighborhood of infinity 
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portions of A and M . Let us set \l/o — on A and X = a = p = 7i = on S,fi = const and 
2 9ab to be a (round) 2-sphere metric with area a A on S. Lewandowski [17] has shown that, 
in the resulting solution, A is an isolated horizon with area a A . Note that in the resulting 
solution \&4 need not vanish in any neighborhood of A, or, indeed, even on A. Hence, in the 
vacuum case, there is an infinite-dimensional space of (local) solutions containing isolated 
horizons. It turns out that, in this setting, there is always a vector field £ a in a neighborhood 
of A with £ a = £ a and ^g a b — 0. However, in general ^Cabcd 0, where C a b c d is the Weyl 
tensor of g a \,. Hence, in general, £ a cannot be a Killing field of g a b even in a neighborhood of 
A. (For details, see |T7[.) The Robinson-Trautman solutions provide interesting examples 



of exact solutions which bring out this point f20fl : a sub-class of these solutions admit an 
isolated horizon but no Killing fields whatsoever. There is radiation in every neighborhood 
of the isolated horizon but, in a natural chart, the metric coefficients and several of their 
radial derivatives evaluated at A are the same as those of the Schwarzschild metric at its 
event horizon. 

A second class of examples can be constructed by starting with Killing horizons and 
'adding radiation'. To be specific, consider one asymptotic region of the Kruskal space-time 
(figure 4) and a Cauchy surface M therein. The idea is to change the initial data on this 
slice in the region r > 3m, say, where m is the Schwarzschild mass of the initial space-time. 



In the Einstein-Maxwell case, one can use the strategy introduced by Cutler and Wald [21 



in their proof of existence of solutions with smooth null infinity. In the vacuum case, one 



can use the more general 'gluing methods' recently introduced by Corvino and Schoen 22 



to show that there exists an infinite-dimensional space of asymptotically flat initial data on 
M which agree with the data for a Schwarzschild space-time for r < 3m but in which the 
evolved space-time admits radiation. Using these methods, one would be able to construct 
'triangular regions' M. bounded by M, A and a partial Cauchy slice M' in the future of M. If 
one takes Ai as the space-time of interest, then A would serve as the isolated horizon at the 
inner boundary. Due to the presence of radiation, M. will not admit any global Killing field. 
However, in a neighborhood of A, the 4- metric will be isometric to that of Schwarzschild 
space-time. Thus, in this case, there will in fact be four Killing fields in a neighborhood of 
the isolated horizon. 

The two constructions discussed above are complementary. The first yields more gen- 
eral isolated horizons but the final result is local. The second would provide space-times 
which extend from the isolated horizon A to infinity but in which there is no radiation in 
a neighborhood of A. We expect there will exist an infinite-dimensional space of solutions 
to the vacuum Einstein equations and Einstein-Maxwell equations which are free from both 
limitations, i.e., which extend to spatial infinity and admit isolated horizons with radiation 
arbitrarily close to them. However, a comprehensive treatment of this issue will be tech- 
nically difficult. Given the current status of global existence and uniqueness results in the 
asymptotically flat contexts, the present limitations are not surprising. Indeed, the situation 
at null infinity is somewhat analogous: while the known techniques have provided several 
interesting partial results, they do not yet allow us to show that there exists a large class of 
solutions to Einstein's vacuum equations which admit complete and smooth past and future 
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null infinities, J?^ and the standard structure at spatial infinity, i°. 



4 Consequences of the Boundary Conditions 

In this section, we will discuss the rich structure given to the horizon by the boundary 
conditions. The discussion is divided into four subsections. The first describes the basic 
geometry of an isolated horizon. The next two subsections examine the restrictions on the 
space-time curvature and Maxwell field which arise from the boundary conditions. The last 
subsection contains a brief summary. 



4.1 Horizon Geometry 

Let us begin by examining the consequences of condition III in the main definition. The 
condition on oa immediately implies 

VJb = -2UJ b (4.1) 

for some 1-form U a on A. Hence, £ a is geodesic, twist-free, shear-free and divergence- free. 
We will denote the acceleration of £ a by k: 

tVJ b = k£ b , (4.2) 

so that k = —2£ a U a . Note that k varies with the rescaling of £ a . 

Actually, the geodesic and the twist-free properties of £ a follow already from condition 
II which requires A to be a null surface with £ a as its null normal. Furthermore, since 
T ab £ a £ b > by condition |V4 and Einstein's equations hold at A (condition IVb), we can 
use the Raychaudhuri equation 

JS$0 W = -\9 2 {l) - a ab a ab + uj ab u ab - R ab t£ b (4.3) 

to conclude the shear o ab must vanish if the expansion 9m vanishes. Thus, the only inde- 
pendent assumption contained in the first of equations ( |3.4j ) in condition III is that l a is 
expansion- free, which captures the idea that A is an isolated horizon. Finally, note that the 
Raychaudhuri equation also implies that 

R ab t£ b = MGT ab £ a £ b = (4.4) 

i.e., that there is no flux of matter energy-momentum across the horizon. 



Let us now consider the properties of the vector field n a . The second equation in (|3.4j ) 
implies 

Van b = (2U a n b + 2/im (a m fe) ), (4.5) 

with U a oc n a . Hence, n a is twist-free, shear-free, has spherically symmetric expansion, 2/x, 
and vanishing Newman-Penrose coefficient 7r = rh a £ b 'V a n b . The twist-free property follows 
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from the very definition of n a , but the other three properties originate in condition III of the 
Definition. 

Next we turn to the intrinsic geometry of the horizon. Using the definitions ( |3.3| ) of the 
co- vectors m a and fh a , one can show there exists a 1-form v a on A such that 

dm = —iv A m and dm = iv A m. (4-6) 

One consequence of these relations is that the Lie derivative along i a of the intrinsic metric 
on A must vanish: ££tg a b — 0. Thus, as mentioned in Section 3, the intrinsic geometry 
of an isolated horizon is 'time-independent'. In particular, the Lie derivative along l a of 
the volume form on the foliation 2-spheres vanishes: \ — 0. Therefore, the areas of 
all the iSa take the same value which we denote a A . Finally, one can show |23|] the scalar 
curvature 2 R of the 2-sphere cross sections is related to the 4- dimensional, Newman- 
Penrose curvature scalars via: 2 R = — 2Re [NE^] + 2$n + R/12. In Section [5], it is shown that 
the quantity on the right side of this equation is constant on Hence, the intrinsic metric 
2m( a ?fife) on Sa is spherically symmetric. In this sense, the horizon geometry is undistorted. 
However, the discussion of Section [T2] shows spherical symmetry will not extend, in general, 
to a neighborhood of A. 



4.2 Form of the Curvature 

Let us begin by exploring the effects of boundary condition |V] on the form of the Ricci 
tensor. Using the Raychaudhuri equation, we have derived (|4.4j). Whence k a in ( |3T7| ) must 
be proportional to £ a . Using the quantity e defined in ( |3.8|) , we then have 



(SttG)- 1 (R ab t - \Ri a + M a ) = T ab £ b = -e£ a , (4.7) 

where R is the scalar curvature. This formula yields a series of results for the Ricci tensor. 
In terms of Newman- Penrose components ( |A.14j ), these read 

R A 

$oo = $oi = $10 = and $ u + - - - = AnGe. (4.8) 

o 2 

The first three results say, by way of the Einstein equation, there is no flux of matter radiation 
falling through the isolated horizon. The fourth result implies the combination $n + -| is 
spherically symmetric. 

We can now explore the consequences of the condition III for the full Riemann curvature. 
Since any SL(2, (D)-bundle over a 3-manifold is trivializable, and since our 4-manifold M. has 
the topology M x R for some 3-manifold M, the connection A can be represented globally 



as a Lie algebra- valued 1-form. Because of ( |3.4[) , in the (t, o)-basis, the pull-back to A of 
that connection must have the form 

A a AB = -{kn a + ivy A o B ^ - fim a o A o B . (4.9) 
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on A. (However, since the spin-dyad is defined only locally on A, the decomposition 
is also local.) The function \x appearing here is the same nowhere vanishing, spherically 
symmetric function introduced in ( |3.4| ) and k and v are defined by ( |4.2| ) and (|4.6|), respec- 
tively. Now, using this expression for the pull-back of the self-dual connection, we can simply 
calculate the pull-back of its curvature to be 



AB 



■[dk An + idv\i^OB) + [(^/i + *5?en)n A m]oAOB, 



(4.10) 



where we have suppressed all space-time indices for simplicity. On the other hand, we 
can also calculate the self-dual part of the Riemann spinor in terms of Newman-Penrose 
components (see Appendix ||). Then, using the compatibility of the self-dual connection 
with the soldering form at the horizon, we get a second expression for the pulled-back 
curvature: 



AB 



+ 



$oo n Am + ^/ n Am — (\l>i — §oi)m A m iaIb 

$i n A m + vj^n Am - {^> 2 ~ $11 - t^J m Am 2l^ a o B ) 

$20^ Am + (^5/ 2 + n Am — (\^ 3 — $21)^ A m OaOb- 



(4.1i; 



Equating this expression with ( |4 . 1 U| ) one arrives at a series of conclusions: 
1. Since there is no laLb term in fl4.10|) , we find 

$0 = #i = 0, 

where we have used (|4.8|). 



2. Equating the L(aOb) terms in the two expressions and using ( 4.12 ) and 



n A dn = 



and 



dv = 2i ( — $11 - — ] m Am. 



(4.12) 

then yields 
(4.13) 



The first expression here says the function k is spherically symmetric. In the second 
expression, the left side is real and the quantity im Am = \ on the right is real as 
well. As a result, the coefficient ^2 — $n — || must be real. However, || is manifestly 
real and, due to its definition (|A.14|) , $ n is also real. Thus, the second expression in 
( f4.13| ) implies the imaginary part of ^2 vanishes. This encodes the property that A is 
non-rotating. 

3. Equating the remaining oaOb terms similarly yields $20 — and \l/3 = $21 as well as 



R 

^ 2 + 12 ~ ^ + 



(4.14) 



Since /1 and n are spherically symmetric, ^2 + t| must also have this property. 
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In the later sections of this paper, we will consider the action and phase space formula- 
tions of systems containing isolated horizons. In this discussion, it will be most useful to have 
a formula giving the relations which arise from the boundary conditions among the funda- 
mental gravitational degrees of freedom in a simple, compact form. Using (|4.11| ), the tetrad 
decomposition ( |A.9| ) of T, AB and the above restrictions on the Newman- Penrose curvature 
components, it is straightforward to demonstrate: 



m 2 - $ n 



24 J 



V 2 



$1; 



A B 
O O L C LD 



^CD 



(4.15) 



In the phase space formalism, only the pull-back to Sa of this formula will be relevant. This 
pull-back takes the simpler form 



F AB = (tf 2 

Note that the essential content of this equation can be seen already in (|4.13| ). 



R 



$n E 

11 24 ' 



AB 



(4.16) 



4.3 Form of the Maxwell Field 

The stress-energy tensor of a Maxwell field is given by 

1 



Tab ~ 47T 



F ac F b c -^ a6 F cd F crf . (4.17) 



This stress-energy tensor satisfies the dominant energy condition and hence, in particular, 
condition |V'a| . Furthermore, one can see from its definition that the trace of T a b is zero. 

To see the restrictions which the boundary conditions place on the Maxwell field, it is 
useful to recast this discussion in terms of spinors as we did in the previous subsection for 
the self-dual curvature. The Maxwell spinor (f) A B is defined in terms of the field strength via 
the relation 

F afe = a AA> a h BB \(t> AB e AIBI + £ab4>A' b 1 ) ■ (4.18) 

It is then straightforward to show that the stress-energy (|4.17| ) can be expressed in terms of 
the Maxwell spinor as 

T a6 =-^/ A V B >^<W (4-19) 

We have already seen in the previous subsections that the general matter field conditions, |VJ 
and the Raychaudhuri equation imply the stress energy tensor must satisfy (|4.7|) . Using the 
spinorial definition ( |3.2j ) of £ a , this restriction on the stress-energy tensor gives two important 
results for the Maxwell field: 

0o := <p AB o A o B = and e = ^ AB i A o A ■ 4>A'B'l A 'o B ' = r- \<f>if ■ (4.20) 

Z7T Z7T 



Here, <e is the quantity e introduced in ( |3.8j ) specialized to the Maxwell field. The second 
equation shows the spherical symmetry of we imposed in (|3.5|) does indeed guarantee 
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the spherical symmetry of <e required by ( |3.8|) . However, the remaining Newman- Penrose 
component of the Maxwell field, 02, is completely unconstrained. In particular, it need not 
be spherically symmetric. 

Since fa is spherically symmetric, we can express it in terms of the electric and magnetic 
charges contained within the horizon. To do so, consider the general form of the Maxwell 
field compatible with ( 4.20 ): 



F = I A (0 2 m + farh) - 2 (im [0!] \ - Re [fa] * 2 e) . (4.21) 

Here, 2 e denotes the volume form on Sa and * 2 e denotes its (space-time) dual. Now, since 
\ is defined with respect to the inward pointing unit space-like normal (see Appendix [A]), 
we have 



and 



4ttQa = f *F = -2f Re [fa] \ (4.22) 



-4ttP a = f F = -2f Im [fa] % (4.23) 



where Qa and Pa denote the electric and magnetic charges contained within Sa- Since fa is 
spherically symmetric, its real and imaginary parts can be pulled outside the integrals and 
we calculate it to be 

fa = — (Qa + ^Pa). (4.24) 

«A 

Here, Qa and Pa are naturally spherically symmetric, but may as yet vary from one Sa to 
another. However, the remaining boundary condition on the Maxwell field, { \ Vc|) , requires 



the Maxwell equations hold at the horizon. The Maxwell equations pulled-back to A applied 
to the field strength (|4.21| ) with fa given by (|4.24j ) imply the charges Qa and Pa must be 



constant over the entire horizon A. It should be noted that this constancy is caused by 
boundary conditions and not by equations of motion in the bulk. As a result, Qa and Pa are 
constant over A in any history satisfying our boundary conditions and not just 'on-shell'. 



4.4 Summary 

As we have seen in this section, the boundary conditions place many restrictions on both the 
gravitational and electro-magnetic degrees of freedom. We will collect the results we have 
found here. These results use not only the boundary conditions, but also the fact that the 
only form of matter we consider is a Maxwell field. 

The Newman-Penrose components of the Maxwell field at the horizon are constrained by 

27T 

0o = O and 0i = — (Qa + ^Pa)- (4.25) 

a A 

The remaining component, 02, is an arbitrary complex function over A. 
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The Newman-Penrose components of the Ricci tensor are 



R = 4A and = 2G^j with i,j = 0,1,2. (4.26) 



The second equation is simply a consequence of the Einstein equation and Q4.19 ). The 
Newman-Penrose components of the Weyl tensor satisfy 

•q/ = \fr 1 = and ^ 3 = $ 2 i- (4.27) 

The component ^4 is an arbitrary complex function over A. The remaining component, ^2, 
is related to the acceleration k of i a and the expansion 2fi of n a via: 

* 2 + - = Sffm + Rfi, (4.28) 

where A is the cosmo logical constant. It follows in particular that ^2 is real and spherically 
symmetric. Finally, these components also satisfy 

du = 2 f* 2 - $n - -] % (4.29) 



6 

where v is a connection on the frame bundle of Sa and 2 e is its volume form. 



5 Surface Gravity and the Zeroth Law 

The zeroth law of black hole mechanics states that the surface gravity k of a stationary black 
hole is constant over its horizon. In subsection |5.1| , we will extend the standard definition of 
k to arbitrary non-rotating isolated horizons using only structure available at the horizon. 
A key test of the usefulness of this definition comes from the zeroth law: Does the structure 
of A enable us to conclude k is constant on A without reference to a static Killing field? In 



Section |572| we will show the answer is in the affirmative. Thus, our more general definition 
of surface gravity is consistent with our notion of isolation of the horizon. Furthermore, we 
will see that the structure of A is rich enough to enable us to express k in terms of the 
parameters r&, Qa and Pa of the isolated horizon. 



5.1 Gauge reduction and surface gravity 

Our set-up suggests we define surface gravity using the acceleration, k, of the vector field £ a . 
However, since the acceleration fails to be invariant under the rescalings of £ a , we need to 
normalize l a appropriately. As mentioned in the Introduction, the usual treatments of black 
hole mechanics in static space-times accomplish this by identifying i a with the restriction to 
the horizon of that static Killing field which is unit at infinity. For a generic isolated horizon, 
there will be no such Killing field and our procedure can only involve structures defined on 
the horizon. Since £ a is null, and its expansion, twist and shear vanish, we cannot hope to 
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normalize it by fixing one of its own geometric characteristics. However, the normalizations 
of l a and n a are intertwined and we can hope to normalize n a by fixing its expansion. The 
normalization of l a will then be fixed. 

To implement this strategy, let us begin by examining the available gauge freedom. The 
correspondence ( |3.2j ) between the fixed spin dyad (i A , o A ) and the preferred direction fields 
[f\n a ] breaks the original SL(2, C) internal gauge group at A to C£/(l). However, as we 
shall see shortly, the residual gauge invariance has a somewhat unusual structure. To see 
the nature of the residual gauge, it is simplest temporarily to consider a fixed soldering form 
and a variable spin dyad. 

The most general transformation of the spin dyad which preserves the correspondence 

(PDis 

(i A , o A ) h- (e e - ie t A , e" e+i V) . (5.1) 

Here, G and 9 are both real functions over A. Under these residual gauge transformations, 
one can show the null tetrad transforms as 



t ^ e" 2e r n a i-> e 2e n a 
m a i-> e 2W m a m a ^ e" 2ie m a . 



(5.2) 



Thus, G accounts for the rescaling freedom in i a and n a which must be eliminated to define 
the surface gravity unambiguously. Q is restricted to be spherically symmetric by (|3~l| ). On 
the other hand, 9 is arbitrary and allows a general transformation on the frame bundle of each 
5*a- Furthermore, the function fi appearing in ( |3.4|) is not gauge invariant, but transforms 
according to 

fi i-> e 2e /i. (5.3) 
Note, however, that the spherical symmetry and nowhere-vanishing property of fi are pre- 



served by this transformation. Finally, the fields k and v introduced in Section [4J] transform 

as 

k i-> e- 2& (k -2J2&9) and v a i-> u a - 2V a 9. (5.4) 

The transformation of k is the usual one for the acceleration of a vector field when that field 
is rescaled, and preserves the spherical symmetry of k. Meanwhile, the transformation of v a 
is that of a U(l) connection. 

We are now ready to fix the normalization of i a . The strategy outlined in the beginning 
of this subsection can be implemented successfully thanks to the following three non-trivial 
facts. First, the expansion of the properly normalized n a in a Reissner-Nordstrom solution is 
'universal': irrespective of the mass, electric and magnetic charges or cosmological constant, 
on the future, outer black hole event horizons in these solutions, 9i n \ = —2/tq where r is 
the radius of the horizon. Motivated by this and the relation 6^ = 2fi (see Q3.6D ) we are led 
to setf|| 

H= (5.5) 



4 To accommodate cosmological horizons, we will have to allow fi to be strictly positive. This issue will 
be discussed in section |[ For now, we focus on black hole horizons and assume /i is strictly negative. 
Modifications required to accommodate positive \x are straightforward. 
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on a generic, non- rotating, Einstein-Maxwell isolated horizon of radius r A (i.e., of area 
a A = 47rr A ). Second, the gauge-freedom (|5.2|) available to us is such that we can always 
achieve the normalization ( [5.5| ) of fi. Third, it is obvious from ( |5.3| ) that this condition 
exhausts the freedom in O completely. In particular, therefore, in a single stroke, this 
procedure fixes the normalization of n a (and £ a ) and gets rid of the awkwardness in the 
gauge freedom. The restricted gauge freedom is now given simply by: 

r ^ t n a ^ n a 

m a ^e 2ie m a m a ^ e~ 2ie m\ ( ' J ' (>) 

where 6 is an arbitrary real function on A; the local gauge group is reduced to U{\). We can 
now return to our usual convention wherein the spin dyad is fixed while the soldering form 



varies. The true residual gauge transformations are then the duals of (|5.1|) , with 6 = 0, 
acting on the soldering form. The effect of these residual transformations on the null tetrad 
and the other fields discussed here remain the same as those in (|5.6| ). 

From now on, we will assume that n a and £ a are normalized via ( f5.5| ), denote the resulting 
acceleration of i a by k and refer to it as the surface gravity of the isolated horizon A. 
By construction, our general definition reduces to the standard one in Reissner-Nordstrom 
solutions. 

To conclude this subsection, let us consider the gauge freedom in Maxwell theory. We 
just saw that a partial gauge fixing of the SL(2, C) freedom in the gravitational sector is 
necessary to define the surface gravity in the absence of a static Killing field. The situation 
with the electric and magnetic potentials is analogous. In conventional treatments [^4| 
these can be defined using the global static Killing field which, however, is unavailable for 
a generic isolated horizon. The idea again is to resolve this problem through a partial 
gauge fixing. Since the only available parameters are the radius r A and the charges Qa, -Pa, 
dimensional considerations suggest the electric potential $ = A a £ a is proportional to Qa/^a 
on the horizon. We fix the proportionality factor using the standard value of $ in Reissner- 
Nordstrom solutions. Thus, in the general case, we partially fix the gauge by requiring 

tK a = $ := (5.7) 

r A 

It turns out that this strategy of gauge-fixing also makes the variational principle well-defined 
for the Maxwell action. 

The situation with the magnetic potential, however, is more subtle. Since there is no 
obvious expression for the magnetic potential in terms of A a , we cannot formulate a definition 
similar to fl5.7|) . Instead, we will appeal to the well-known duality symmetry of the Maxwell 
field. Thus, in the remainder of the paper, we set Pa — in the main discussion. The results 
for isolated horizons with both electric and magnetic charge will follow from those including 
only electric charge by a duality rotation. 
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5.2 Zeroth law 



We already know from ( 4.13|) that k is spherically symmetric. Therefore, to establish the 
zeroth law, it only remains to show that Jzf^K also vanishes. Recall first from ( |4.15| ) that on 
A, F AB is given by: 



^9 - $11 - — I 6^5 



$11 A B L C L D 



^CD 



(5.8) 



Let us now consider the Bianchi identity DF = 0. Transvecting it with laOb, we obtain 



&t $2 -$11 



R 
24 



0. 



(5.9) 



In the Einstein-Maxwell system, $ n = 2G|0i| 2 and <f>\ is constant on A (see ( 4.24Q ). Simi- 
larly, since the stress-energy tensor is trace-free, R = 4A is also a constant. Hence it follows 
that =Sf£$2 — 0. Finally, ( |4.28| ) and our gauge condition ( |5\5| ) immediately imply: 



K 



T A 

$2 + — 

3 



(5.10) 



Hence, we conclude J^k = 0, as desired. This establishes the zeroth law of the mechanics 
of isolated horizons in the Einstein-Maxwell theory. 

We will now obtain an explicit expression for k in terms of the parameters of this isolated 
horizon. The key fact is that the first Chern number of the pull-back to 5a of the connection 
v introduced in ( |4.6| ) is two. This can be seen in a number of ways, but it is essentially 
equivalent to the Gauss-Bonet theorem for a 2-sphere because v can be identified with a 
50(2) connection in the frame bundle of 5a- Using this property, we have 



(f idv = — (f 2 

2iri Js A 2tt Js A 



#2 + 



R_ 

12 



$n + 



R 



(5.11) 



where we have used ( |4.13j ) in the last equality here. Now, we have seen in ( f4.14| ) that the 
first term in the brackets here is spherically symmetric, and in (|4.8| ) that the second term is 
proportional to the quantity e which has been restricted to be spherically symmetric. As a 



result, the entire integrand on the right side of |5.11| can be pulled through the integral. The 
remaining integral simply gives the area a A of 5a- Thus, 



R ~ 27r 

$2 - $11 " — = 

24 a A 



(5.12) 



Let us now specialize to the Einstein-Maxwell case. Then, R = 4A and $n = (G/2r^)(Q' A + 
P A ). Therefore, using ( p,10| ) we can express surface gravity in terms of ta, Qa and Pa: 



It, 



(5.13) 
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We conclude with a few remarks. 



1. The final expression (|5.13|) for k is formally identical to the expression for the surface 
gravity of a Reissner-Nordstrom black hole in terms of its radius and charge. This may 
be surprising since we have not restricted ourselves to static situations. However, if it is 
possible to express k in terms of the parameters r&, and A alone, this agreement 
must hold if the general expression is to reduce to the standard one on event horizons of 
static black holes. In our treatment, the agreement can be technically traced back to our 
general strategy for fixing the normalization of l a . 

2. It may also be surprising that, although we do not have a static Killing field at our 
disposal, it was possible to define k unambiguously and it turned out to satisfy the zeroth 
law. Furthermore, we could express k in terms of the parameters of the isolated horizon, 
irrespective of the details of gravitational and electro-magnetic radiation outside the horizon. 
This was possible because of two facts. First, the boundary conditions could successfully 
extract just that structure from static black holes which is relevant for these thermodynamical 
considerations. Second, at its core, the zeroth law is really local to the horizon; it does not 
know, nor care, about the space-time geometry away from the horizon. Physically, this meets 
one's expectation that the degrees of freedom of a black hole in equilibrium should 'decouple' 
from the excitations present elsewhere in space-time. 

3. Our expression ( [4.14] ) of surface gravity in terms of Weyl and scalar curvature is universal, 
i.e., holds independent of the particular matter sources so long as they satisfy the mild energy 
condition (|3.7|) . Furthermore, in all these cases, k is spherically symmetric and the Bianchi 
identity ensures (|5.9| ). The restriction to Maxwell fields as the only source has been used 
here only to show that -^q^h ($ii + f — f) = ^eT ab £ a n b = Jzf^e vanishes on A. Thus, in 
the present setting of non-rotating isolated horizons, the zeroth law would hold for more 
general matter provided its stress energy tensor satisfies this last restriction. This condition 
is satisfied, for example, by dilatonic matter |7], [16]. 



4. In the main definition, we assumed A has topology S 2 x R. If S 2 is replaced by a compact 
2-manifold of higher genus, the results of Section |] and the proof of constancy of k on A 
remain unaffected. However, in obtaining the explicit expression ( |5.13[ ) of surface gravity, 



we used the Gauss-Bonnet theorem. Hence this expression is not universal but depends on 
the genus of 5a ■ 



6 Action and Hamiltonian 

As pointed out in the Introduction, to arrive at an appropriate generalization of the the 
first law of black hole mechanics, we first need to define the mass of an isolated horizon. 
The idea is to arrive at this definition through Hamiltonian considerations. In Section |6.1| , 
we introduce an action principle which yields the correct equations of motion despite the 



presence of the internal boundary A. In Section |6.2j, we pass to the Hamiltonian theory 



by performing a Legendre transform and in 6.3 we discuss Hamilton's equations of motion. 
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We will find that, due to the form of the boundary conditions, there are subtle differences 
between the Lagrangian and the Hamiltonian frameworks because the latter allows more 
general variations than the former. 

As in previous sections, in the main discussion we assume the space-time under consid- 
eration is asymptotically flat with vanishing magnetic charge and discuss at the end the 
modifications required to incorporate non-zero A and Pa- 

6.1 Action 

Recall from Section that, in absence of internal boundaries, the action for Einstein-Maxwell 
theory is given by: 

S(a, A, A) = — ^ / Tr[S A F] + ^ f F A *F + — ^— / Tr [£ A A] (6.1) 

07TG J M 07T J M 07rG JCoo 

In the presence of internal boundaries, however, this action need not be functionally differ- 
entiable. This is the case with our present boundary conditions at A. In 0, the required 
modifications were discussed for the case when all histories (a, A, A) under consideration 
induce a fixed area a A and electro- magnetic charges Qa (and Pa) on A and a differentiable 
action was obtained by adding to S a surface term at A. (For a general discussion of surface 
terms in the action, see P5| .) 

The strategy of fixing the parameters of the isolated horizon was appropriate in |7j because 
the aim of that analysis was to provide the classical framework for entropy calculations of 
horizons with specific values of their parameters. In this paper, on the other hand, we 
need a Hamiltonian framework which is sufficiently general for the discussion of the first 
law, in which one must allow the parameters to change. Therefore, we now need to allow 
histories with all possible values of parameters. Note, however, that a key consequence of the 
boundary conditions is that the area a A and charge Qa of the isolated horizon are constant 
in time. Therefore, the values of these parameters are fixed in any one history (although 
they may vary from one history to another). Now, since all the fields are kept fixed on the 
initial and final surfaces in the variational principle and the values of our dynamical fields 
on either of these surfaces determine a A and Qa for that history, one is only allowed to use 
those variations for which 5a A and 5Qa vanish identically. This fact simplifies the task of 
finding an appropriate action considerably: For example, as far as the action principle is 
concerned, one can continue to use the action used in 0. 

There is however, a further subtlety which will lead us to use a more convenient boundary 
term in the action. Because of the nature of the variational principle discussed above, we 
are free to add any function of the horizon parameters a A and Qa without affecting the 
Lagrangian equations of motion. In the framework considered in J/J, this just corresponds to 
the freedom of adding a constant (with appropriate physical dimension) to the Lagrangian. 
However, since the full class of histories now under consideration allows arbitrary areas a A 
and charges Qa, the freedom is now more significant: it corresponds to changing the La- 
grangian by a function of parameters a A , Qa- While the variational procedure used to derive 
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Figure 5: Region Ai of space-time considered in the variational principle is bounded by 
two partial Cauchy surfaces M\ and M 2 . They intersect the isolated horizon A in preferred 
2-spheres Si and S 2 and extend to spatial infinity i°. 



the Lagrangian equations of motion is completely insensitive to this freedom, the Hamilto- 
nians resulting from these Lagrangians will clearly be different. Can all these Hamiltonians 
yield consistent equations of motion? The answer is in the negative. The reason lies in 
subtle differences between the Lagrangian and Hamiltonian variations. In the Hamiltonian 
framework, the phase-space is based on a fixed space-like 3-surface M. Since the values of 
a A and Qa can vary from one history to another, they can also vary from one point of the 
phase space to another. In obtaining Hamilton's equations, 5H = Q(S,Xh), we must now 
allow phase space tangent vectors S which can change the values of parameters a A , Qa- Con- 
sequently, with boundary conditions such as ours, the burden on the Hamiltonian is greater 
than that on the Lagrangian. It turns out that these additional demands on the Hamiltonian 
suffice to eliminate the apparent functional freedom in its expression. More precisely, the re- 
quirement that Hamilton's equations of motion be consistent for all variations 5 in the phase 
space, including the ones for which 5a A and 5Qa do not vanish, determine the Hamiltonian 
completely. (There is no freedom to add a constant because, with only Newton's constant 
G and the speed of light c at our disposal, there is no constant function on the phase space 
with dimensions of energy.) One can then work backwards and single out the expression of 
the action, which, upon Legendre transform, yields the correct Hamiltonian. We will follow 
this strategy. 

Fix a region of space-time whose inner boundary is an undistorted, non-rotating isolated 
horizon, as depicted in figure 5. This region M. is bounded in the past and future by 
space-like hypersurfaces Mi and M 2 respectively which intersect the horizon A in preferred 
2-sphere cross-sections Si and S% and extend to spatial infinity i°. Since the space-time M. 
is asymptotically flat at spatial infinity, the fields obey the standard falloff conditions at 
%° . The interior boundary A is a non-rotating isolated horizon which satisfies the boundary 
conditions discussed in section 3. It turns out that, to obtain a well-defined action principle, 
we need to impose an additional condition at A: 
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for any (preferred) 2-sphere cross-section Sa of A, where \ is the natural alternating tensor 
on Sa (see Appendix A. 2). Note that this restriction is very mild since it only asks that the 
spherically symmetric part (with respect to 2 e) — or, equivalently, the zero mode — of A - £ 
be real on A. Then, the required action is given by: 



S =lTr! t Tr[SAF] + -/ F A *F + f Tr[E A A] + — I "(r A * 2 ) % (6.3) 
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where A e is the volume form on A compatible with our normalization for £ a (see Appendix 



A). In Section |6T3| , we find the resulting Hamiltonian does yield a consistent set of equations. 
That discussion will also show that the term added in the passage from ( |6.1| ) to ( |6.3j ) is 
uniquely determined by the consistency requirement. 

Note that ( |6.3| ) does not have the Chern-Simons boundary term at A used in p]. However, 
if one restricts oneself to histories with a fixed value of a a as in [0], ( |6.3|) is completely 
equivalent to the action used there. (The mild restriction ( |6.2|) was not discussed in but 
is needed also in the action principle used there.) In particular, as we will see in Section 



672] , the symplectic structure obtained from the present action ( |6.3D again has a boundary 
term at A. If one works with a fixed a A , this term reduces to the Chern-Simons symplectic 
structure of @. For non-perturbative quantization || §, it is this symplectic structure that 
plays the important role; the form of the action is not directly relevant. 

In this paper we have chosen to work with ( |6.3|) because it is more convenient for the 
Hamiltonian framework with variable a A , needed in the discussion of the first law. Further- 
more, this form of the action appears to extend naturally to isolated horizons with distortion 
and rotation and also may be better suited for quantization in these more general contexts. 



6.2 Hamiltonian Framework 

To pass to the Hamiltonian framework, we need to perform the Legendre transform. As in 
Section [2.2| , we begin by introducing the necessary structure on the 4-manifold J\4. First, 
foliate J\A by partial Cauchy surfaces M with the following properties: i) the 3-manifolds M 
intersect the horizon at the preferred 2-spheres Sa such that the unit time-like normal r a to 
them coincides with the vector (£ a + n a )/\^2 at Sa', and, ii) they extend to spatial infinity 
and intersect in 2-spheres S^. Next, fix a smooth time-like vector field t a , transverse to 
the leaves M and a function t such that: i) t a V a t = 1 on Ai; ii) t a tends to the unit time 
translation orthogonal to M at spatial infinity; and, iii) t a tends to l a on the horizon A. 
(The restriction on t a that it be orthogonal to the foliation at infinity has been made only 
for simplicity and can be removed easily by suitably modifying the discussion of the physical 
interpretation of surface terms in the Hamiltonian.) Finally, we will restrict ourselves to 
physically interesting situations in which M are partial Cauchy surfaces for the space-time 
region M. under consideration and adapt our orientations to the case in which the projection 
of n a into M is a radial vector which points away from the region M.. (See figure 6 and the 
discussion that follows in Section |6.3j ). 
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The 3+1 decomposition of the space-time fields can now be performed exactly as in 
equations (|2.6p and ( |2.7|) . Once again, the phase space consists of quadruples (A a AB , 
S ab j4s ,A a ,E a f ) ) on the 3-manifold M satisfying appropriate boundary conditions. The condi- 
tions at infinity are the same as before, namely ( [2.8|) . However, there are additional boundary 
conditions at the horizon. First, because of ( f4.9p , the form of the gravitational connection 
A is restricted at 

A AB = -iv 1^ + —m o A o B . (6.4) 
r A 

Next, the curvature of A is restricted by (|4.16 ) and (|5.11| ) to satisfy F = — 27rS/a A , and 



the electro-magnetic field F is restricted by ( 4.24 ). Finally, the requirement that the action 



principle be well-defined imposes the mild restriction (|6.2j ) at A. 

The Legendre transform is straightforward using the procedure outlined in Section ^72 
The only new element is the treatment of boundary terms at the horizon which requires the 
use of the boundary conditions listed above. In order to state the final result, we have to 
introduce some notation. The space of our connections v on A has the structure of an affine 
space. Let us choose any one of these connections v, satisfying v ■ i = as the 'origin'. 
(For example, v can be the 'static magnetic monopole potential' on every S A -) Then using 
boundary conditions, it is easy to show that any other connection v can be expressed as: 

v = v + r] + dip (6.5) 

where the 1-form rj on A satisfies: i ■ rj = 0, i ■ drj = and d*r] = where the dual is taken 
under the metric on each Sa- Note that there is the obvious gauge freedom ip i— > ip + const 
in the choice of the function ip. Let us eliminate it by requiring 

ip 2 e = (6.6) 

5 A 



on any Sa- Then the Legendre transform of the action S of ( |5.3|) is given by: 

S(a,A,A)= fdt ^— [ J? t AAE--p- <f ^\ ~ [ dtH t , (6.7) 
J [8nG JM(t) 8tvG Js A (t) J 

with 

H t = H t -f (I^^-^AX (6.8) 

where H t is defined in ( |2.9| ). Thus, the Hamiltonian has the familiar form: As in ( |2.9| ), the 
bulk term is a volume integral of constraints weighted by Lagrange multipliers determined 
by t a and the surface term at infinity gives jf; a p ADM = —£j ADM . However, now there is now 
a surface term at the horizon as well. 

In order to bring out the similarity and differences in the two surface terms, let us express 
the term at infinity using the Weyl curvature. Assuming the field equations hold near infinity 
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and with the shift N set to zero at Soo because of our current assumptions on the asymptotic 
value of t a , we have 



H t = I constraints + lim i ( ^ 2 ) \ - J> ( * 2 - ) \ (6.9) 

Jm r o-"^Js ro \47rG / Js A \4ttG 27rr A J 

Thus, only the 'Coulombic parts' of the two curvatures enter the expressions of the two 
surface terms. However, while the surface term at infinity depends only on the gravitational 
curvature, the term at the horizon depends also on the Maxwell curvature. 
The symplectic structure also acquires a term at the horizon. 



fiOMa) = V($i,$2) - ^ / M<5 2 ( 2 e) - 8 2 1) 6, ( 2 e)] (6.10) 
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where f2 is defined in ( |2.9| ). (As mentioned in Section p.l| , if we restrict ourselves to the 
phase space corresponding to a fixed value of a&, the surface term in (|6.10j ) reduces to 
the Chern-Simons symplectic structure for the connection A on Sa)- The presence of the 
surface term in the symplectic structure is rather unusual. For instance, although there 
is a boundary term in the action at infinity, the symplectic structure does not acquire a 
corresponding boundary term. Also note that we have not added new 'surface degrees of 
freedom' at the horizon (in contrast with, e.g., |27], |28|). Indeed, our phase space consists 
only of the standard bulk fields on M which normally arise in Einstein-Maxwell theory (see 
Section |2]2|) and whose values on 5a are determined by their values in the bulk by continuity. 
If anything, the boundary conditions restrict the degrees of freedom on A by relating fields 
which are independent in the absence of boundaries. The symplectic structure on the space 
of these bulk fields simply acquires an extra surface term. In the classical theory, while this 
term is essential for consistency of the framework, it does not play a special role. In the 
description of the quantum geometry of the horizon and the entropy calculation [|| §], on 
the other hand, this term turns out to be crucial. 

6.3 Hamilton's Equations 

Hamilton's equations are 

SH t = n(6, X H ), (6.11) 

where Xh is the Hamiltonian vector field associated with the given time evolution vector 
field t a and 5 is an arbitrary variation of the fields. Unlike in the discussion of the action, all 
fields appearing in the Hamiltonian are defined only on the space-like surface M. Hamilton's 
equations describe the time evolution of these fields. In particular, there is no a priori reason 
to expect the area or charge of the isolated horizon to be constant in time. The constancy 
of the area and charge must arise, if at all, as equations of motion of the theory. As we 
already noted, since the linearized fields 5 in ( |6.11|) can have Sa A ^ and 5Qa ^ 0, there are 



'more' Hamilton's equations than what one would have naively thought from the Lagrangian 
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perspective. The question is whether the additional equations ensure the area and charge 
are conserved and if the full set of equations is self-consistent. 

Let us summarize the consequences of Hamilton's equations for the Hamiltonian and 
symplectic structure introduced in the last subsection. The bulk equations of motion give 
the standard Einstein-Maxwell equations expressed in terms of connection variables. As 
usual, the variation of the term at infinity in the Hamiltonian cancels the surface term at 
infinity arising from the variation of the bulk term. Also, the equations of motion preserve 
the boundary conditions at infinity. 

Thus, it only remains to examine the horizon terms. Using the relation ( |5.1U| ) among 
\I/ 2 , k and r A , and equating the horizon terms on the two sides of Hamilton's equations, we 
obtain: 

2 e = and ip = vl. (6.12) 

The first of these equations in particular guarantees that the horizon area does not change 
under time-evolution. The second equation follows from (|6.5|) which defines ip and is thus a 
consistency condition. Note also that the restriction (|6.6|) on ip is preserved in time because 

of (|3). 

Finally, it is natural to ask whether Hamilton's equations imply Qa — 0. The answer is 
in the affirmative. However, this result is a consequence of a bulk equation of motion which 
guarantees § s E = where 5* is any closed two surface. If we take S = Sa, the obvious 
consequence is 

Qa = 0. (6.13) 

To summarize, there exists a unique consistent Hamiltonian formulation in the presence 
of inner boundaries which are isolated horizons. The symplectic structure is given by ( |6.10| ), 
and the Hamiltonian by ( |6.9| ). The bulk equations of motion are the standard 3+1 versions 
of the Einstein-Maxwell equations. There are, however, additional equations on the horizon 
2-sphere Sa which guarantee that r A is a constant of motion. 

For simplicity, in the main discussion we restricted ourselves to zero magnetic charge 
and cosmological constant, only one asymptotic region and only one inner boundary. How- 
ever, these restrictions can be easily removed. If there is more than one asymptotic region 
and/or isolated horizon inner boundary, one need only include surface terms for each of these 
boundary 2-spheres. The incorporation of a non-zero magnetic charge and cosmological con- 
stant has a slightly more significant effect. As in Section 2^, the presence of a cosmological 
constant changes the boundary conditions and the surface terms at infinity. The symplectic 
structure is unchanged. But, as discussed in Section 7TT , the surface term at the horizon in 
the expression of the Hamiltonian acquires additional terms involving P A and A. 



7 Physics of the Hamiltonian 



In this section, we will examine the expression (|6.9| ) of the Hamiltonian in some detail and 
extract physical information from it. In Section 771 , we will argue that the surface term 
at Sa should be identified with the mass of the isolated horizon. In Section |7.2|, we will 
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show that the numerical value of the Hamiltonian in any static solution vanishes identically 
so that the mass of the isolated horizon in such a space-time reduces to the ADM mass at 
infinity. Finally, in Section [7.3| , we will argue that, in any solution to the field equations 
which is asymptotically flat at null and spatial infinity and asymptotically Schwarzschild at 
future time-like infinity, the Hamiltonian ( p.9\) equals the total energy radiated away through 
future null infinity, J? + . In these space-times, the mass of the isolated horizon then equals 
the future limit of the Bondi mass, exactly as one would intuitively expect. We should 
emphasize, however, that the argument rests on assumptions on the asymptotic behavior of 
various fields (particularly near i + ) and we do not prove the existence of solutions to field 
equations with this behavior. Therefore, the discussion of Section |7]3| has a different status 
from the rest of the paper. Its primary purpose is to strengthen our intuition about the 
Hamiltonian and the mass of the isolated horizon. 



7.1 Isolated Horizon Mass 

For any physical system, energy can be identified with the numerical, on-shell value of the 
generator of the appropriate time translation. In Minkowskian field theories, for example, it 
is the generator of motions on phase space which correspond to space-time diffeomorphisms 
along a constant time-like vector field. Consider, as a second example, the theory of gravita- 
tional and electro-magnetic radiation in general relativity. In this case, one can construct a 
phase space of radiative modes at null infinity and the total radiated energy is the numerical 
value of the Hamiltonian generating a time translation in the BMS group at null infinity 



30|j . Finally, in the physics of fields which are asymptotically flat at spatial infinity, 
the ADM energy arises as the on-shell, numerical value of the Hamiltonian generating an 
asymptotic time-translation. If the space-time under consideration admits several asymp- 
totic regions (as, for example, in the Kruskal picture) then the energy in any one asymptotic 
region is given by the generator of a diffeomorphism which is an asymptotic time-translation 
in the region under consideration and asymptotically identity in all other regions. 

These considerations suggest we define the energy associated with a given isolated horizon 
A to be the numerical, on-shell value of the generator of a diffeomorphism which is a time 
translation at A and asymptotically identity. To obtain an expression for this energy, let 
us examine the expression ( |6.9|) of the Hamiltonian H t . The bulk term vanishes on shell 
and the term at infinity does not contribute if the vector field t a vanishes at spatial infinity. 
Thus, the required expression is given simply by the surface term at Furthermore, since 
the vector field t a tends to i a on A and, by construction, i a defines the 'rest-frame' of the 
isolated horizon, this energy can be identified with the mass Ma of A. Thus, the Hamiltonian 
considerations lead us to set 



M A = - 

/5 A 



2. 
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where we have now allowed for a non-zero cosmo logical constant A and magnetic charge Pa- 
For purposes of the first law, it will be useful to rewrite this expression by eliminating 
the curvatures ^2 and <\>\ in favor of surface gravity k and electro-magnetic scalar potential 
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$. Using flOl ) and (|5TT0|) , when P A = 0, we have: 



M A = -^-kaa + ^Qa (7.2) 
47rG 



Note that the expression is formally identical with the familiar Smarr formula [[31] for the 
mass of a Reissner-Nordstrom black hole. One also knows directly (i.e., without making an 
appeal to black hole uniqueness theorems) that the ADM mass of any static black hole in 
the Einstein-Maxwell theory is given by (|7.2|) ||24j| . Thus, as with our definition of surface 
gravity k, although Ma is defined using only the structure at the isolated horizon A, it agrees 
with the standard definition of black hole mass for static solutions. The reason behind this 
agreement will become clear in the next sub-section. However, in general (non-static) space- 
times, due to the presence of radiation, the ADM mass at infinity is quite distinct from the 
isolated horizon mass Ma- When constraints are satisfied, we have 

H t = M A - E ADM , (7.3) 

and we will see in Section I7T31 that the numerical value of —H f can be identified with the total 



energy radiated through future null infinity. Finally, note that Ma has a specific contribution 
from the Maxwell field. We will see that this contribution is rather subtle but quite crucial 
to adequately handle charged processes in the first law. As far as we are aware, none of the 
general, quasi-local expressions of mass contain this precise contribution from the Maxwell 
field. Thus, in the charged case, it appears that Ma does not agree with any of the proposed 
quasi-local mass expressions. 

A natural question is whether Ma is positive. Let us first consider the case with zero 
cosmological constant. Then, by fixing the value of the charge Qa and minimizing Ma with 
respect to r A , one finds Ma is not only positive, but also bounded below: M\ > GQ\. 
At the minimum, i.e. when M\ = GQ\, the surface gravity k vanishes. However, unlike 
Ma, k can be negative. This structure is familiar from Reissner-Nordstrom solutions, where 
the same inequality holds, k vanishes at extremality, is positive on the outer horizon and 
negative on the inner. However, it was not obvious that this entire structure would remain 
intact on general isolated horizons. 

Let us now consider the case when the cosmological constant A is non-zero. If A is nega- 
tive, Ma is again positive and M A > GQ\. In this case, is natural to impose asymptotically 
anti-de Sitter boundary conditions, whence one only expects 'black-hole type' horizons. If A 
is positive, one also has cosmological horizons and the situation becomes more involved. The 
resulting complications are illustrated by the Schwarzschild-de Sitter space-time (see figure 
6). The Hamiltonian framework is physically useful only in those situations in which M 
is a partial Cauchy surface for the space-time region hA under consideration. The isolated 
horizons in this case are future boundaries of space-time such as (a) and (b) in the figure. 
Surface gravity as well as mass are positive on the black hole horizon (a). The case of the 
cosmological horizon requires a reconsideration of the sign conventions we previously adopted 
(see footnote 4). Specifically, in the construction of the Hamiltonian framework of Section 
|6.1|, we chose our orientations by assuming the projection of n a into M is 'outward pointing' 
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Figure 6: Schwarzschild-de Sitter space-time. Surface (a) is black hole horizon, (b) and (d) 
are cosmological horizons, and (c) is a white hole horizon. The isolated horizon boundary 
conditions are satisfied on all four. However, the Hamiltonian framework is tailored to future 
horizons of types (a) and (b) which are boundaries of space-time regions admitting partial 
Cauchy surfaces M. The expansion of n a (and hence /x) is negative on (a) and positive on 
(b). Surface gravity k and mass Ma are positive on both (a) and (b). 



at Sa relative to M. With this choice, the expansion of n a (and hence the Newman-Penrose 
coefficient fi) is negative on the black hole horizon (a), but positive on the cosmological 
horizon (b). Since we assumed, again for definiteness, that \i is negative in Sections || and 
|6], certain trivial modifications are needed to accommodate cosmological future horizons of 
the type (b). With these changes, the surface gravity and mass are again positive on (b). 

To summarize, it is future horizons of type (a) and (b) that are of physical interest 
in our Hamiltonian framework. For them, the surface gravity and mass are positive in 
Schwarzschild-de Sitter space-time and we expect the situation to be similar for general 
isolated horizons with positive cosmological constant. More detailed considerations suggest 
that the interpretation of a A as entropy and k as temperature are applicable only to such 
horizons. 

7.2 Static solutions 

The phase space now under consideration admits a 2-parameter family of static solutions, 
labeled by M and Q — the Reissner-Nordstrom solutions. Let us begin by evaluating the 
Hamiltonian H t (of fl6.9| )) on these solutions using for t a the static Killing field. Then, the 
volume integrals will vanish since the constraints are satisfied and only contributions from 
Soo and Sa will remain. The term at infinity equals the negative of the ADM mass while, as 
noted above, the horizon term is given by + $Qa- Now, it is well known from the Smarr 
formula that on a Reissner-Nordstrom space-time, the value of the ADM mass is precisely 



33 



+ $Qa- Therefore the value of the Hamiltonian H t when evaluated on static space-times 
is zero. (The same reasoning extends to the case when A and Pa are non-zero.) 

This feature is not accidental; there is a general argument from symplectic geometry 
which 'explains' this vanishing of H t . We will conclude this sub-section by presenting the 
argument. 

In symplectic geometry, Hamilton's equations are 5H = Q(5, Xh), where 5 is an arbitrary 
variation and Xh is the Hamiltonian vector field. A stationary solution (such as a Reissner- 
Nordstrom solution in the sector of Einstein-Maxwell theory now under consideration) is one 
at which the Hamiltonian vector field either vanishes or generates pure gauge evolution. In 
either case, the symplectic structure evaluated on Xh and any other vector field 5 vanishes. 
Therefore, at these points of the phase space, 5H = for any variation 5. In particular 
5H = for variations relating two nearby stationary solutions. Let us suppose the phase 
space is such that the space of stationary solutions is connected (an assumption satisfied 
by the Reissner-Nordstrom family in our case). Then, the Hamiltonian must take some 
constant 'preferred value' on all stationary solutions. Now, let us suppose there is no natural 
energy scale in the theory. (This assumption is satisfied in our case because Ma and Qa are 
not fixed on our phase space and because one cannot construct a quantity with dimensions 
of mass from G and c alone.) Then, this 'preferred value' must be zero. 

7.3 Hamiltonian equals Radiative Energy 

We now present a result which provides an intuitive interpretation of the Hamiltonian H t 
and a further motivation for our definition of the isolated horizon mass Ma- More precisely, 
using suitable regularity assumptions, we will show that, when the equations of motion 
(with A = 0) are satisfied and t a is adapted to the natural rest frames at the horizon, the 
numerical value of H t equals t a P^ a,d = —E™ d , where E™ d is flux of energy radiated across 
J^ + . However, as explained in the beginning of Section [F|, because we will need a number of 
new assumptions, the considerations of this sub-section are not as self-contained as those of 
the rest of the paper. 

Let us assume that the underlying space-time M. is of the type indicated in figure (l.a). 
That is, we assume i) the space-time is asymptotically flat at future null infinity ,J? + and 
asymptotically Schwarzschild at time-like infinity i + ; ii) the Bondi news tensor on ^ + tends 
to zero as one approaches i° and i + ; iii) the isolated horizon A extends to i + ; and, iv) the 
boundary of A4 consists of A, U i + U i° and a partial Cauchy surface M. 

Fix a conformal completion {M. 1 g ab ) 1 of (M.,g a b) which has J? + as its (future) null 
boundary. Appendix C summarizes the structure available at Let us begin by recalling 
that part of its structure which we will need in this subsection. The conformal factor Q 
vanishes at and h a := V a f2 is the null normal to J^ + . The conformally rescaled metric 
g a b induces a degenerate metric q ab at which satisfies q a bV a = on if and only 
if the tangent vector V a to J? + is proportional to n a . Although is degenerate, we can 
define its 'inverse' q ab via q a bQ bc Qcd = (lad- This 'inverse' is unique up to additions of terms 
of the form n^ a V b ^ for some vector field V a tangent to J? + . The phase space of radiative 
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modes of the gravitational and electro-magnetic fields at can be coordinatizedQ by pairs 
of fields (7 a fe, A a ) defined intrinsically on J^ + . The fields 7 a b code the gravitational degrees 
of freedom; they are symmetric, transverse (i.e., 7 a &n 6 = 0) and trace- free (i.e., r ) a bQ ab = 0). 
These properties imply that 7 has precisely two independent components which represent 
the two radiative modes of the gravitational field. The Maxwell degrees of freedom are coded 
in the 1-form fields A a on J^ + , satisfying A a n a = 0, with A tending to zero at i°. Again, A a 
has two independent components which capture the two radiative modes of electro-magnetic 
field. The symplectic structure can be written as: 



fi md (5l, S 2 ) : = / [Sllab ^n^lcd) ~ ^lab %n{hlc d )\ 

+ ±- [ q ab [SiA a ^n(S 2 A b ) - 6 2 A a J4(5i A b )] y e 



(7.4) 



(For details, see Appendix C.l and p{| |3~0 



The asymptotic symmetry group at is the BMS group |32] which admits a preferred 
four- dimensional Abelian sub-group of translations. Let us suppose that the conformal factor 
is chosen such that q ab is the unit 2-sphere metric. Then, h a is a BMS time-translation. 
Diffeomorphisms generated by n a induce motions on the radiative phase-space. As one 
might expect, they preserve fi rad and the corresponding Hamiltonian is given by ||30|| : 



H™ d = / N ab N cd q ac q M ^ - — / (F ac F b c + *F ac *F 6 c )n a n 6 J e (7.5) 

where N ab = —IS^fCiab is the Bondi News tensor at the point in the radiative phase space 
labeled by 7^. Thus, 5H rad = f2 rad (<5, Xn) for any tangent vector 5 to the radiative phase 
space. Using the asymptotic form of the space-time metric in suitable coordinates, Bondi 
and Sachs |32| identified the right side of (|7.5|) as 



P rad . t = -E r&a (7.6) 

where P rad is the 4-momentum radiated across J !+ and t represents the BMS time translation 
defined by n a . Thus, _£7 ad is the flux of energy across carried by gravitational and electro- 
magnetic waves. (Again, the negative sign arises in ( |7.6| ) because our signature is — h ++.) 
The Hamiltonian formulation at null infinity 0] provides a general conceptual setting in 



support of this interpretation. 

We wish to relate these structures on the radiative phase space with those on the canonical 
phase space introduced in Section |6]. Fix a point on the constraint surface of the canonical 
phase space and evolve it using field equations. Consider tangent vectors satisfying linearized 

5 More precisely, the radiative phase space consists of certain equivalence classes of connections on J> + . 
Thus it has a natural affine space structure. In introducing this coordinatization, an equivalence class of 
'trivial' connections has been chosen as the origin. The fact that the phase space is an affine — rather than 
a vector — space has some subtle but important consequences. These will be ignored here as they do not 
affect the issues now under discussion. For details, see [E9| [30[. 
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constraints at this point and evolve them using linearized field equations. In appendix C, 
assuming the background and linearized solutions satisfy certain falloff conditions, we find 



fi(<$i, 6 2 ) = fi rad (<5r, % )■ ( 7 - 7 ) 

where f2 rad is evaluated at the point in the radiative phase space defined by the background 
solution and S Tad is the tangent vector in the radiative phase space defined by the linearized 
solution 5. The idea is to let 5\ be arbitrary, choose for 62 the Hamiltonian vector field 
defined by a time translation, and use ( |7.7|) to relate the canonical Hamiltonian (|6.9| ) to the 
radiative Hamiltonian ( [7.5|) . 

Let us choose a vector field t a on M. such that: i)t a = £ a on A; ii) t a is unit at spatial 
infinity and defines an asymptotic time-translation; and, iii) t a is a BMS time translation at 
^ + and the conformal factor is so chosen that t a = n a at J^ + . As in Section || leaves M 
of our foliation will be assumed to be asymptotically orthogonal to t a . Then, from Section 
I we have 5H t = Sl(6,X t ), with X t = (^f t A,J^E; 3f t A,£tE). From the above discussion of 
the radiative phase space, we have: 5 rad H™ d = fi rad (5 rad , X[ ad ) where X[ ad = (Jzf A7 ; J&A). 
Therefore, using the equality ( |7.7| ) of the two symplectic structures, we conclude: 

SH t = 5 rad H™ d = -5 rad E™ d (7.8) 

for all linearized solutions 5 satisfying the asymptotic conditions. Let us assume such so- 
lutions span the tangent space at every point in the portion of the phase space under con- 
sideration. Then, we conclude that H t and if t rad differ by a constant. To fix the value 
of that constant, let us examine the Reissner-Nordstrom space-times. We already saw in 



Section |7.2| that in these space-times H t vanishes. Since they are static, the Bondi News 
tensor and the electro-magnetic radiation vanish on J? + . Hence, if[ ad also vanishes in a 
Reissner-Nordstrom solution. Thus, the value of the constant is zero and 

Ht = -E™ d . (7.9) 

As one might have intuitively expected, the canonical Hamiltonian equals the time compo- 
nent of the flux of the 4- momentum that is radiated across J? + . 

Recall from ( |7.3| ) that, when the constraints are satisfied, the value of the canonical 
Hamiltonian H t is equal to the isolated horizon mass minus the ADM energy E ADM in 
the rest frame of the isolated horizon. Therefore, it now follows that Ma = E ADM — E™ d . 
It is well-known that the difference of the ADM energy and the flux of energy through 
denoted E Bondl (i + ), is the future limit of the Bondi energy (in the rest frame defined by t a ) 
||33| . Hence, we conclude: 

M A = E Bondi (i + ) (7.10) 

Thus, Ma can be thought of as the mass remaining in the space-time after all radiation has 
escaped to infinity, or, equivalently, the mass of the black hole with its static hair. This 
simple interpretation provides additional support for our definition of the horizon mass. 

As emphasized earlier, the discussion of this subsection is based on a number of tech- 
nical assumptions (which are stated in Appendix C.2). We will conclude with a summary 
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of their physical content. Apart from asymptotic flatness at spatial and null infinity, the 
key assumptions involve the structure of i + . We assume there is only one bound state in 
asymptotic future, represented by the isolated horizon. This reflects the expectation that, in 
the Einstein-Maxwell theory, there would be no gravitational or electro-magnetic radiation 
hovering around the horizon at late times. Multiple black hole solutions which reach equi- 
librium asymptotically are excluded, as their structure at i + would be quite different from 
that of the Schwarzschild space-time. If the black holes do not reach asymptotic equilibrium 
but accelerate away from each other, the structure at i + may be similar to (or even simpler 
than) that in the Schwarzschild space-time. An example is provided by the C-metric, where, 
if the parameters are adjusted suitably, the structure at i + as well as that at i° is regular 
as in Minkowski space-time. However, the accelerating black holes pierce J" + which is now 
singular. Therefore, as it stands, our analysis is not applicable to this case either. Although 
our analysis could conceivably be generalized to cover these two types of situation, its cur- 
rent form is primarily applicable to the situation depicted in figure 1(a) in which a single 
gravitational collapse occurs. 



8 First Law 

Since we now have well-defined notions of surface gravity k, electric potential $ and horizon 



mass M A , we are ready to examine the question of whether the first law holds. In section [8]1 
we will consider the equilibrium version of the law in which the horizon observables of two 
nearby space-times are compared. This version is closer in spirit to the treatment of the first 
law of thermodynamics in which one compares the values of macroscopic, thermodynamic 
quantities associated with two nearby equilibrium configurations, without reference to the 



process which causes the transition between them. In section |8.2| , we will consider the physical 
process version of the first law in which one explicitly considers the process responsible for 
the transition. This version will bring out certain subtleties. 

8.1 Equilibrium version 

Denote by XH. the (infinite-dimensional) space of space-times admitting (one or more) iso- 
lated horizons. In this section, we will be concerned only with the structure near isolated 
horizons. In particular, we will not have to refer at all to the boundary conditions at infinity 
or to the precise nature of the matter outside the isolated horizon. We will simply assume the 
surface gravity k, the potential $ and the mass M A of the isolated horizon are determined 
by its intrinsic parameters r A and Qa via ( |5.13| ), ( |5.7| ) and (|7T2|): 



MA = ^(l + ^%-Ar A 2 ) a A =4vrr A 

... _ 1 I 1 G< 2a , 2 \ *_ ^ 



2 ^A V r A J r A 
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and allow general matter fields in the exterior. (Recall from Section |^ that the notion of 
isolated horizons is not tied down to Einstein-Maxwell theory.) This viewpoint is similar 
to that normally adopted for the ADM 4-momentum and angular momentum defined at 
spatial infinity. These quantities are first derived from Hamiltonian considerations adapted 
to specific matter sources (e.g., Klein Gordon, Maxwell, Dirac and Yang-Mills fields) but then 
used also for total 4-momentum and angular momentum for more general forms of matter 
(e.g. fluids) for which a satisfactory initial value formulation and Hamiltonian framework 
may not exist. Thus, one often uses the expressions of 4-momentum and angular momentum 
at infinity without specifying the precise matter content, assuming only that the stress- 
energy tensor satisfies physically reasonable conditions and falls off appropriately. In the 
same spirit, we now assume that Ma, k and $ of an isolated horizon A of radius r A and 
charge Qa are given by ( |8.1| ), irrespective of the matter content outside, so long as that 
matter does not endow the horizon with additional intrinsic parameters (such as a dilatonic 
charge or a new U(l) charge). 

Given a space-time (A4, g ab ) in TH and a tangent vector 5 at this point, we can consider 
a smooth curve in TH passing through this point with 5 as a tangent vector there and 
examine how M A and k associated with the isolated horizon A in the background change. 
Straightforward algebra yields: 

SM A = -^K5a A + $6Q A . (8.2) 
87rG 

This equation tells us the relation between infinitesimal changes in the mass, area, and charge 
of two nearby, non-rotating isolated horizons. It is our generalization to isolated horizons of 
the equilibrium version of the first law of black hole mechanics. 
We will conclude this sub-section with a few remarks. 

1. The above calculation leading to the first law is trivial. The non-trivial part of the 
analysis was to arrive at expressions (|8.1| ) of k and Ma in absence of a static Killing field. 
Again, although our boundary conditions allow the presence of radiation arbitrarily close to 
the horizon, they successfully extract the structure from event horizons of static black holes 
that is relevant for thermodynamic considerations. As with the zeroth law, the veracity of 
the first law can be taken as additional support for our definitions of k and Ma- 

2. The laws of black hole mechanics were first derived by Bardeen, Carter and Hawking 
|2|. |3|. They considered stationary black holes possibly surrounded by a perfect fluid in a 
circular flow and arrived at the first law by comparing two nearby stationary solutions. For 
purposes of comparison, it is more convenient to use an extension of that work to more 



general matter fields discussed by Heusler [24]. In the non-rotating case, their main results 



can be summarized as follows. Identities governing the Komar integral of the static Killing 
vector imply the mass Madm measured at spatial infinity is given by 

Madm = ~r^a h + f (2T ab - Tg ab ) K a dS b (8.3) 
47rG Jm 
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where is the area of the horizon, M is a partial Cauchy surface from the horizon to spatial 
infinity, T a ^ is the stress-energy of the matter field, and K a is the static Killing field. Under 
variations 5 from one such stationary black hole solution to another, it was then shown that 

5M ADM = « SA+ 1 j Qab{5 y K _ J_ f S*(G.K) , (8.4) 

87TG 167T JM 87T JM 

where *if is the 3-form dual to K a and *(G.K), the 3-form dual to G ab Kb. In the Einstein- 
Maxwell case, it turns out that the two volume integrals in ( |8.4j ) collapse to a single term 
at the horizon which is precisely <MQa> where $ = A.K. Thus, in this case, the first law 
relates the change in the ADM mass to changes in quantities at the horizon. The final form 
is the same as ( |8.2| ) (with Ma replaced by Madm)- Similarly the final expression of Madm 
is just the right side of ( [7.2|) . 

However, there are some important differences from our approach. First, in the above 
derivation, one restricts oneself to static solutions of field equations and makes a heavy use of 
the Komar integral associated with the Killing vector K a . Second, the permissible variations 
6 are only those which relate nearby static solutions. Third, the term §5Qa at the horizon 



arose from the volume terms in (O); that its value depends only on the fields evaluated at 



the horizon is not fundamental to this derivation. The whole calculation of |24| is based on 
an interplay between infinity and the horizon which is possible only because the assumption 
of staticity makes the problem elliptic and field variations 'rigid'. The black hole is not 
studied as an isolated, separate entity; the quantities defined at the horizon are tightly tied 
to the exterior fields. Indeed, there is no useful analog of the mass Ma, associated with 
the horizon. Perhaps the closest analog is the Komar integral evaluated at the horizon, 
sometimes denoted Mh @- However, in the charged case, this integral does not include 
the Maxwell contribution $<5a and cannot therefore be used directly in the first law. Our 
derivation, by contrast, makes no reference to a Killing vector and allows radiation fields 
outside the horizon. All our considerations are local to the isolated horizon; in the variation, 
we did not have to refer to bulk fields on M at all. Finally, the mass which appears in our 
first law is the isolated horizon mass, Ma, and not the ADM mass Madm- 

3. A treatment of the first law based on Hamiltonian considerations was given by Wald 



and collaborators |34| , p5| , p6fl . The final result of this work is more general than that of 
Bardeen- Carter-Hawking type of analyses. The background space-time is again a stationary 
black hole, possibly with matter fields in the exterior. However, the perturbations 6 are no 
longer required to be stationary; they can relate the background stationary solution to any 
nearby solution. Since our approach is also based on the Hamiltonian framework, the two 
treatments share a number of common features. 

However, there are also a number of important differences, both in methodology and 
final results. While boundary conditions play a key role in our analysis, their analogs are 



not specified in references |34], |35|, |36|. Consequently, the issue of differentiability of action 
is not discussed. In particular, while our action S of ( |6.3| ) cannot be written as a pure 
bulk term, in p3, BR |36| there is no surface term in the action either at infinity or at the 



horizon. Consequently, the Hamiltonian contains only bulk terms and there is no analog of 
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our horizon mass Ma- When restricted to non- rotating black holes, the first law of |34|. |35|. [36 



has the same form as ( |S.2|) . However, as in the Bardeen-Carter-Hawking approach, our Ma 
is replaced by the ADM mass Madm and the background space-time is assumed to be a 
stationary black hole solution. 



4. As indicated in Section |3.2j , the space TTi of space-times admitting isolated horizons in 
Einstein-Maxwell theory is infinite-dimensional. The space S of static solutions is a finite- 
dimensional subspace of TH. In the Bardeen-Carter-Hawking type approach, the first law 
holds only at points of S and for tangent vectors 5 to S. In the Wald approach, it holds 
again only at points of S but the variations 5 need not be tangential to S. In the approach 
developed in this paper, the first law holds at all points of 17i and for all tangent vectors 
to XJi. Thus, the generalization involved is very significant. However, since our boundary 
conditions imply the intrinsic geometry of 5a is spherically symmetric, distorted black holes 
are excluded from our analysis. By contrast, the other two approaches can handle static, 
distorted black hole solutions. In Einstein-Maxwell theory, there are no such solutions. 
However, if we allow charged fluid sources, such solutions presumably exist. Therefore, in 
the general context, our framework misses out certain situations which are encompassed by 
the other two approaches. It would be interesting to generalize our framework to overcome 
this limitation. 



8.2 Physical process version 

Let us now consider the situation depicted in figure 1(b). We are given a space-time with a 
non- rotating, isolated horizon Ai (with parameters t\ and Q\). Suppose a small amount of 
matter falls in to the horizon and after a brief dynamical period the horizon settles down to a 
new equilibrium configuration A 2 (with parameters r 2 and Q<i). The question is: How do the 
observables associated with the horizons change in this physical process? The difference from 
the situation considered in the previous subsection is that one is now considering a physical 
process occurring in a single space-time rather than comparing two nearby space-times. 

It is completely straightforward to analyze the process in our framework since the two 
masses and surface gravities are determined by their intrinsic parameters via (|8.1|) . The 
actual algebraic calculation is the same as in the last sub-section; only the physical meaning 
of the variation S is different. Hence, we again find the changes in mass, area and charge are 
governed by ( |S.2| ). 

It is instructive to analyze this relation in terms of properties of the matter which fell 
across the horizon. Using the Raychaudhuri equation, and keeping only first order terms in 
variations, it is straightforward to show |36| 

K 5a A = 87rG 5E Rux = 8tiG [ *(5T-£), (8.5) 

J 

where is the portion of the horizon (between Ai and A 2 ) crossed by the matter and 
*(5T ■ £) is the 3-form dual to 5T a b£ b . Let us first suppose Q 2 = Qi, i.e., the charge of the 
horizon did not change in this physical process. Then, comparing (|8.5| ) with (|8.2j ), we arrive 
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at a simple physical picture: the change in the mass of the horizon is equal to the total energy 
flux E anx across the horizon. However, it is interesting to note that, if Q2 7^ Qi, there is an 
extra contribution, $5Qa, to What is the origin of this term? It arises due to 'book- 

keeping' in the following sense. As we observed in Section [7|, Ma contains not just the 'raw 
energy of the content of the horizon', but also the energy of the electro-magnetic hair outside 
the horizon. (Recall that in static solutions, Ma equals the ADM mass and, more generally, 
it equals the future limit of the Bondi mass, both of which include the contribution from 
energy in the Coulombic electro- magnetic field outside the horizon.) Before the physical 
process began, the charge 5Q& = Q2 — Qi is outside the horizon and the energy in its 
Coulomb field does not contribute to Mi. At the end of the process, however, the black hole 
charge changes by 5Qa and the energy in the corresponding Coulombic field does contribute 
to M2. This accounts for the term $<5Qa in the expression of £Ma = M2 — Mi. Thus, 
the physical process version of the first law is subtle. The first order change in the mass of 
the horizon has a two-fold origin: a contribution due to flux of energy across horizon and 
another contribution from book-keeping of the energy in the Coulombic hair of the horizon. 

What is the situation in the standard framework, where one uses Madm i n place of Ma? 
To our knowledge, the physical process version of the first law has been discussed only in 
the uncharged case ^6f . One assumes that the background space-time is globally static and 
considers a (non-static) matter perturbation which falls across the horizon. The ADM mass 
of the unperturbed space-time is taken to be Mi and the ADM mass of the background 
plus perturbation is taken to be M 2 . Then, using the reasoning given above, one arrives at 
the first law SM = (nSa) /87rG and interprets SM as the change in the mass of the black 
hole due to the energy flux across the horizon. However, in the charged case, if is not 
equal to zero, it seems difficult to account for the term QSQa which also contributes to 
SM without bringing in Ma-Q It is interesting to note that, in this respect, there is a key 
difference between the angular momentum work term Q5J and the electro-magnetic term 
$ 5Q/\: While the angular momentum contribution is coded easily in the flux of the stress- 
energy across the horizon, the electro- magnetic contribution is not. This is why, unlike the 
electro-magnetic work term, the angular momentum work term can be easily incorporated 
in the physical process version in the standard approach [RH . 



9 Discussion 

Let us begin with a summary of the main ideas and results. 

In Section |3] we introduced the notion of a non-rotating isolated horizon A. While one 
needs access to the entire space-time to locate an event horizon, isolated horizons can be 
located quasi-locally. Event horizons of static black holes in Einstein-Maxwell theory do 

6 At first sight, it may appear that one should be able to account for this term using Maxwell's equations 
to simplify the contribution to J« * (ST ■ I) that comes from the first order change in the stress energy of the 
Maxwell field. However, a closer examination shows that the integral of the Maxwell contribution vanishes 
identically. 
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qualify as isolated horizons. However, the definition does not require the presence of a Killing 
field even in a neighborhood of A. Rather, physically motivated, geometric conditions are 
imposed on the null normal i a to A and on an associated inward pointing null vector n a at 
A. These conditions imply the Lie derivative along l a of the intrinsic (degenerate) metric 
of A vanishes which in turn implies the area of an isolated horizon is constant in time. In 
this sense, the horizon itself is isolated or 'in equilibrium'. However, the space-time may 
well admit electro-magnetic and/or gravitational radiation. The quasi- local nature of the 
definition of A and the possibility of the presence of radiation suggest the space of solutions to 
the Einstein-Maxwell equations admitting isolated horizons would be m/imte-dimensional, 
in striking contrast to the space of static black holes which is only three-dimensional. Recent 
mathematical results by a number of workers |[7|, [18|, 19, |22| show this expectation is indeed 
correct. 

While the conditions used in the definition seem mild, they lead to a surprisingly rich 
structure. In particular, the intrinsic metric of A, the shear, twist and expansion of i a and 
n a and several of the Newman-Penrose gravitational and electro-magnetic curvature scalars 
at A have the same functional dependence on the radius r A and charges Qa,Pa as in the 
Reissner-Nordstrom solutions. This rich structure enables one to fix naturally the scaling of 
the null normal i a and leads to an unambiguous definition of surface gravity, k. Furthermore, 
using only the structure available at A, one can show that k is constant on A; the zeroth 
law is thus extended from static black holes to isolated horizons. 

To formulate the first law, we need a notion of the mass Ma of the isolated horizon. Since 
we allow for the presence of radiation outside A, we cannot use the ADM mass Madm as Ma, 
nor do we have a static Killing field to perform a Komar integral at A. Fortunately, we can 
use the Hamiltonian framework. Although the presence of the internal boundary introduces 
several subtleties, a satisfactory Hamiltonian framework can be constructed. When the 
constraints are satisfied, the Hamiltonian turns out to be a sum of two surface terms, one 
at infinity and one at A. As usual, the term at infinity yields the ADM energy and we 
define Ma to be the surface term at A. This definition is supported by several independent 
considerations. In particular, under suitable conditions, Ma turns out to be the future limit 
of the Bondi mass. Having expressions for both k and Ma at our disposal, we ask if the first 
law holds. The answer is in the affirmative for both the 'equilibrium state' and the 'physical 
process' versions. This provide a significant generalization of the first law of mechanics of 
static black holes in the Einstein-Maxwell theory. Furthermore, in the charged case, this 
analysis brings out some subtleties associated with the 'physical process' version. However, 
since our framework focusses on isolated horizons and small perturbations thereof, it does 
not shed new light on the second law of black hole mechanics which refers to fully dynamical 
situations. 

These underlying ideas overlap with those introduced in references H |7j . In || , Hayward 
introduced, and very effectively used, the notion of 'trapping horizons'. Our isolated horizons 
are a special case of trapping horizons, the most important restriction being our assumption 
that the expansion of the horizon is zero. This assumption is essential to capture the notion 
that the horizon is in equilibrium, which underlies the zeroth and the first law. Furthermore, 
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our method of defining the surface gravity k and the mass Ma of isolated horizons differ 
from those used by Hayward for trapping horizons and consequently our treatment of the two 
laws is also different. (To our knowledge, in the context of trapping horizons, a satisfactory 
definition of surface gravity is available only for spherically symmetric space-times.) However, 
the notion of isolated horizon is clearly inadequate for the treatment of dynamical situations 
which are considered, for example, in the second law and it is these situations that provide 
a primary motivation in the analysis of trapping horizons. 

The relation between the ideas discussed in this paper and those introduced in is 
closer. Both papers deal with isolated horizons. However, while the focus of reference is 
on constructing a Hamiltonian framework suitable for quantization and entropy calculations, 
the focus of the present paper is on the mechanics of isolated horizons. The two overlap in 
their constructions of Hamiltonian frameworks. However, as explained in Section [], reference 
only considers isolated horizons with fixed parameters r A ,Q A ,P A and therefore ignores 
several subtleties which are critical to our present treatment of the first law. Reciprocally, in 

0, significant effort went into the construction of a Hamiltonian framework in terms of real 
variables which is necessary for quantization but not for the laws of mechanics. Finally, in 
Sections |3| and £| and in Appendices A and B, we took the opportunity to present the necessary 
background material from a perspective which is different from but complementary to that 
adopted in 0. 

We will conclude by indicating a few avenues to extend the present work. 

1. Let us begin with the non- rotating case. Although we did not explicitly require the intrin- 
sic metric of an isolated horizon to be spherically symmetric, our assumptions on properties 
of the null vector fields i a and n a at the horizon led us to this conclusion. The discussion 



of Section |2.2| shows the assumptions are not overly restrictive: the class of space-times sat- 
isfying them is infinite-dimensional.[] Typically, these space-times will admit radiation and 
will not be spherically symmetric in the bulk. Nonetheless, it is of interest to weaken our 
assumptions to allow space-times with 'distorted' horizons on which the intrinsic geometry 
will not be spherical. For simplicity, consider the case in which there is no matter in a small 
neighborhood of A. Then, we would expect only to have to weaken the conditions on n a and 
allow fi to be non-spherical. The structure at event horizons of static, distorted black-holes 
has been recently examined by Fairhurst and Krishnan and their analysis confirms this hy- 
pothesis. The extension of the framework presented here to incorporate distortion should be 
fairly straight forward. 

2. Inclusion of rotation would provide an even more interesting extension. Again, conditions 

7 This may seem surprising at first since most of the current intuition comes from static black holes and, 
in the static context, generalizations of Reissner-Nordstrom solutions naturally lead to distorted horizons. 
However, this is because static problems are governed by elliptic equations and generic perturbations in the 
exterior then force the horizon itself to be distorted. Radiative space-times provide generalizations in quite 
a different direction. Now, the equations are hyperbolic and, as the 'gluing methods' of Corvino and Schoen 
p2f show, the geometry can be spherical even in a neighborhood of the horizon without being spherical 
everywhere. More generally, as the Robinson-Trautman solutions illustrate |2(J , the rotational Killing fields 
may not extend even to a neighborhood of the horizon. 
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on £ a would remain unchanged. Only the conditions on n a and the sphericity requirement on 
the component T^l -^ of stress-energy tensor at A will have to be weakened. In particular 
the Newman-Penrose spin-coefficient tt can no longer be zero since it is a potential for the 
imaginary part of ^ which carries the angular momentum information. Work is already in 
progress on this generalization. 



3. In the stationary context, using Hamiltonian methods and Noether charges, Wald |35 



has extended the notion of entropy and discussed the first law in a wide variety of gravita- 
tional theories, possibly coupled to bosonic fields, in any space-time dimension. It would be 
very interesting to extend the present framework for isolated horizons in a similar fashion. 
As a first step, one would recast the framework in terms of tetrads e°[ and the associated 
real, Lorentz connections A a j J . The extension of the resulting (Einstein-matter) action and 
Hamiltonian framework to higher dimensions should then be straight-forward. The first step 
is easy to carry out since the tetrads can be easily obtained from soldering forms and the 
Lorentz connection is just the real part of our self-dual connection. Thus, the 4-dimensional 
tetrad action is, in effect, just the real part of (6.3) and the corresponding Hamiltonian is 
just the real part of ( |6.9| ). Hence, it should be rather easy to extend the present results to 
higher- dimensional general relativity, possibly coupled to matter. Furthermore, since the ba- 
sic variables are tetrads rather than metrics, it should be straightforward to allow fermionic 
matter as well. Incorporating general gravitational theories, on the other hand, could be 
highly non-trivial. 
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A Conventions 

In this paper, capital primed and unprimed indices represent SL(2, C) spinors fields. As 
usual, spinors with only unprimed indices are also interpreted as SU(2) spinors in the phase 
space framework. The spinor conventions are largely those of |[23|| , but with minor modifi- 



cations to replace the H signature of p3| with the — h ++ signature used here. We 

describe these modifications here. 

A.l Metric and Null Tetrad 

The metric is given in terms of the soldering form by 

AA' i\ 1 \ 

9ab = cr a a bAA ,. (A.l) 
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Because of our choice of signature, however, the soldering form must be taken to be anti- 
Hermitian: o a AA = —cr a AA . 

A spin dyad (t A , o A ), satisfying l a oa = 1, defines a null tetrad as follows: 

t = ia a AA , o A o A ' m a = ia a AA , o A l A ' 

(A 2) 

n a = ia a AA , l a I A ' m a = ia a AA , i A o A ' . 

This tetrad obeys the usual inner product conventions in the — h ++ signature: the only 
non- vanishing inner products are £ a n a = — 1 and m a m a = 1. (Note that the definitions of 
m and fh differ from those used in J7|. This change was necessary to make the values of 
our Newman-Penrose curvature components the same as those found in the literature even 
though our signature is — h ++.) 

A. 2 Volume Forms and Orientations 

The volume form on space-time is defined by its spinor expression, which is the same as that 



used in 23 



4 AA 1 BB 1 CC' DD f r * 1 / a o\ 

e abcd = a a a b a c °d l ie AB e CD e A'C' e B'D>+ c - c -\- (A-3J 

This volume form can be expressed in terms of the null tetrad as 

^abcd = 24 ^ l\ a n b m c ™>d\ ■ ( A - 4 ) 

The conventions for inducing volume forms on sub-manifolds of space-time are designed 



to be compatible with those used in [12] and with the usual orientation conventions used in 
Stokes' theorem on Riemannian manifolds. Specifically, this means that a volume form is 
induced on a space-like sub-manifold of space-time by contracting its future-directed, unit 
normal with the last index of 4 e. Then, within a space-like hypersurface, a volume form is 
induced on a two-dimensional sub-manifold by contracting its outward-bound, unit normal 
with the first index of the volume form on the hypersurface. All other orientation conventions 
can be determined from these two. In particular, 



(A.5) 

r in \bc= ^m [b rh c] (A.6) 
-3 \ ab n c] = ~6i n [a m b rh c] . (A.7) 

Here, 3 e denotes the induced volume form on a space-like hypersurface, \ denotes the volume 
form on one of the Sa, and A e denotes the preferred alternating tensor on the null surface 
A. Meanwhile, f denotes the future-directed future normal to the space-like hypersurface 
and ff n denotes the unit radial vector directed inward at the horizon. Note that the inward 
normal is appropriate because Sa is the inner boundary of the space-like surface M. 
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A.3 Self-Dual Basis 

The soldering form defines a basis of self-dual 2-forms on space-time via 

V AB « AA' B r. A' (A B) ( \ q\ 

Kb = 2(J [a a b] A'= 2(y a K °b A'- ( A - 8 ) 

Using the spin dyad, these self-dual 2-forms can be expressed as 

Y> ab AB = A£ [a m b] i A i B + 4 (m [a m b] - £ [a n b] ) i {A o B) - 4n [a m b] o A o B . (A.9) 
One can check that these are, in fact, self-dual in that *Y, AB — iH AB . 

A. 4 Newman-Penrose Components 

We define the Riemann curvature tensor to be 

R abc d k d = 2W [a W b] k c . (A.10) 
The space-time curvature spinors are defined by decomposing the Riemann tensor as 

D „ AA 1 „ BB' „ CC „ DD' 

Kbcd = a a °b a c °d 

X ^A'B'tCD&ABC'D' + ^A'B'^C'D' ^ ' ABCD ~ —R^{A{C^D)B) + C.C.j , (A. 11) 

leading to the expression for the Ricci tensor: 

Kb = °a AA> °b BB ' {~2<5>ABA>B> + ~£W A , B ,} . (A. 12) 

Since the Ricci tensor is real, $>aba'b' is Hermitian. 

In terms of the curvature spinors, we define the Newman-Penrose components of the 
Weyl tensor C a bcd by 

= ^ABCDO A o B o c o D =C abcd £ a m b £ c m d 
*i = ^ABCDO A o B o c t D =C abcd £ a m b £ c n d 

^2 = ^abcdo a o b l c l d =C abcd £ a m b m c n d (A.13) 

^3 = ^ABCDO A L B l C L D = C abcd t U b m C U d 

^4 = ^ABCDi A i B iPi D =C ahcd m a n b m c n d . 

Note that these definitions are the same as those found in the literature 0, |37| and, despite 
the difference in signature, the functions \l/ n take their usual values in specific space-times. 
Similarly, the expressions for the Newman-Penrose components of the Ricci tensor read 



$00 


= $ABA>B>0 A B A 


o B ' 


2 ^ob^" ^ 


$22 


,A B T A' 
— ®ABA>B>t i i 


I B ' 


= \R ah n a n b 




$01 


^ A B-A 

= <Paba>b'0 o o 


I 3 ' 


= \R ab l a m b 


$21 


^ A B-A 
= ^ABA'B'i t O 


T B ' 


= \R a bm a n b 




$02 


i A R A' 
= <$>ABA'B>0 A O a L A 


l B ' 


= \R ab m a m b 


$20 


= ®ABA>B>L A L B d A 


o B ' 


= \R ab m a m b 


(A. 


$10 


= <&ABA>B'0 A L B d A ' 


o B ' 


= \R ah l a rh b 


$12 


= §ABA>B>0 A L B l A 


l B ' 


= \R ab m a n b 





$n = <S>ABA>Bio A L B d A 'l B ' = \R ah (l a n b + m a m b ). 
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As before, these are the standard spinorial definitions for the However, their expressions 
in terms of the Ricci tensor differ from those of [23]| by a minus sign. This difference occurs 



because our Ricci tensor is the negative of the one used in |23| . 

A similar decomposition can be performed on the electro-magnetic field. The Maxwell 
spinor is defined by expressing the field strength as 

F ab = o AA ' a b BB '((j) AB e AIB > + cWa'bO (A. 15) 

Then, the Newman- Penrose components of the Maxwell field are defined by 

00 = <Pabo a o b = -tm b ¥ ab 

01 = <pABt A o B = -\{tn b - m a rh b )F ab = \m a m b (¥ - i*F) ab (A.16) 

02 = <t>ABi A i B = n a m b F ab . 

As with the gravitational field, the values of these functions will be the same as those found 
in the literature. 

Finally, the Newman- Penrose spin-coefficients used in this paper are given by: 

{i = m a rh b V a n b , \ = rh a rh b V a n bl n = tm b V a n b) 

a = -m a m b V a £ b , p = -m a m b V a £ b , e + e = -tn b V a l b . 

Note that, as is common in the black- hole literature, we denote the surface gravity by k (so 
that our k equals (e + e) in the Newman- Penrose notation.) We never need to refer to the 
Newman-Penrose spin coefficient k. 



B Newman-Penrose Components and Self-Dual Cur- 
vature 

The purpose of this appendix is to establish the relation between the self-dual SL(2, (D) 
curvature used in [F2 and the Newman-Penrose curvature components described in | |23| |. 

In any putative space-time where only the Gauss law (|2.2|) is solved, the self-dual curva- 
ture F AB is equal to the self-dual portion of the Riemann curvature defined by 

+ D AB 1 CC _ DD' d AA' BB 1 - m i \ 

n ab — 2 a a a b n CC'DD' e A'B> \ D - 1 ) 

If we now substitute for the Riemann spinor in this expression using ( |A.ll ), one can rearrange 
the terms to yield 

+p AB If A'B' faAB l v CD <\,AB ^ v AB /n 9 \ 

K ab - ~2^ab 9 A'B' ~2 ah W CD~ T^ab ■ V=>.L) 
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It is now a long, but straightforward, process to break this formula into spinor components 
using (|A.9|) . Then, using the definitions QA.13 ) and ( |A.14j ) of the Newman- Penrose compo- 
nents, one can express the result as 



AB 



+ $21)^ An- fy A £ Am- $22^ A m + 

(R \ 
\& 2 + n Am - (\& 3 - $21)^ A rh 

R \ 

q, 2 + $ n _ j £ /\ n _ Am- $ 12 £ A m + 



+ 



$io«- A m + ^>in Am — ^\J/ 2 — $11 — m Am 

(R \ 
^2 + — J £ A m - $02^ A m + 



2i( A o B ) 



$oo^ A m + ^q 72 Am — (^ — $oi) m A m 



:b.3) 



This expresses the self-dual curvature in terms of the Newman- Penrose components in a spin 
dyad satisfying l a oa = 1- The null tetrad here is defined, of course, by the same dyad. 



C Symplectic Structure at Null Infinity 



In section | we used the Legendre transform to introduce a symplectic structure ( |6.10|) on 
the canonical phase space. On the other hand, there is also a natural symplectic structure on 
the space of radiative modes of the Einstein-Maxwell system, defined intrinsically at (future) 
null infinity ,J? + . In this appendix, using field equations, we will show the two symplectic 
structures are equal in an appropriate sense, provided the fields under consideration have 
suitable asymptotic behavior. 

Throughout this discussion, we will restrict ourselves to the region M. of figure 1(a) 
which has as its future boundary and which admits partial Cauchy surfaces M which 
extend from the isolated horizon A to spatial infinity %° . As in Section [7.3| , we will set the 
cosmological constant A to zero. 



C.l Phase Space of Radiative Modes at 

Fix an asymptotically flat space-time (M.,g a b) and consider its Penrose completion (Ai, g a b)- 
As usual, gab = Q 2 9ab is the conformally rescaled metric and is the future null boundary 
of M. where the conformal factor Q vanishes. All fields appearing with a 'hat' will refer to 
the geometry defined by the conformally rescaled metric g a b which is smooth at J^ + . 

Let us begin by recalling the 'universal structure' at null infinity of asymptotically flat 
space-times. First, is topologically S 2 x R. Second, the conformally rescaled metric 
naturally defines an intrinsic, degenerate metric q a b = gab and a null normal field n a = V a fi 
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on J^ + . We shall assume the conformal factor Q is so chosen as to make divergence 
free in the sense that V a n a = on J^ + . By construction, n a defines the unique degenerate 
direction of the intrinsic metric q ab : n a q ab = 0. Therefore, the 'inverse metric' is unique only 
up to the addition of a term of the type h^ a v b ^ where v b is an arbitrary vector field on J^ + . 
(Irrespective of the choice of v a , we have q a b<t c <lcd — Qad)- Finally, the volume 3-form on 
can be defined as 

where 4 e abcd is the volume 4- form defined by the rescaled metric g ab . These structures are 
universal in the sense that they are common to all asymptotically flat space-times; they do 
not carry any information about, e.g., the radiation field which can vary from one space-time 
to another. 

Note however that there remains a conformal freedom in the rescaled metric g ab . If is 
an allowable conformal factor which makes divergence-free, so is by Q' = ujQ, where iv 
is nowhere vanishing on and Jzf^cj = at J zr+ . Under this transformation, the conformal 
metric is rescaled as g ab — > uj 2 g ab . As a consequence, the pairs, (g a b,^ a ) and (u 2 q ab , uj~ 1 h a ) 
are to be regarded as (conformally) equivalent at J^ + . 

We can now turn to the dynamical structures and introduce the radiative modes. Note 
first that the derivative operator V defined by the metric g ab on M. naturally induces a 
derivative operator D defined intrinsically on J* + via the pull-back 

DaKjb = y a K b , (C.2) 

where K_ a is an arbitrary co-vector field defined intrinsically at and K b is any extension 
of Kj, to M.. Since V is metric compatible, it follows that D a q bc = and D a h b = 0. The 
radiative modes of the gravitational field in general relativity are fully encoded in connections 
D on .y + satisfying the above conditions. Recall, however, that there is a residual conformal 
freedom at J^ + . As a consequence, one is led to introduce an equivalence relation between 
connections. The phase space of radiative modes consists of these equivalence classes. It 
thus has the structure of an affine space. The difference between any two connections in 
different equivalence classes can be encoded in a tensor field 7 a b which satisfies 

la b q ab = l ab n b = 7a6 = 7(«6)- (C3) 

Therefore, by fixing a point in the phase space as the 'origin', we can label any other point by 
the corresponding tensor field ^ ab . It is easy to see that ^ ab has two independent components 
which represent the two physical degrees of freedom of gravitational radiation. 

The radiative degrees of freedom of the electro-magnetic field can also be described by 
fields intrinsic to J^ + . It turns out that F ab is completely characterized by the unique 
connection A a at satisfying 

A a n a = and lim A a = fC.4) 

where u is the affine parameter along n a . The connection A a satisfying the above conditions 
has two independent components. These represent the two radiative degrees of freedom of 
the Maxwell field. 
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Thus, the phase space of radiative modes at consists of pairs (7^, A tt ) satisfying the 



(C.5) 



conditions ( |C.3|) and ( |C.4! ) respectively. The symplectic structure on this phase space is 

rr d (^ ad , 5™ d ) : = -J— [ q ac q bd [hia b ^n(S 2lcd ) - 5 2lab J^fccd)] 
+ _L f q^AaJZ^At,) - 5 2 A a ^(5 1 A b )} * r e. 

07T J,J?+ 



For further details, see [29, 30| . 



C.2 Equality of Symplectic Structures 

Let us now return to the canonical phase space of Section || Fix a point on the constraint 
hypersuface and consider tangent vectors which satisfy the linearized constraints. Evolve 
these fields using the appropriate field equations. Then, assuming the resulting 4-geometry 
and the linearized fields thereon satisfy appropriate falloff conditions, they would provide 
a point in the radiative phase space at J? + and tangent vectors at that point. Using this 
correspondence, we will now show the canonical symplectic structure (16.10] ) associated with 
a partial Cauchy surface M (of figure 1(a)) equals the radiative symplectic structure (|C.5| ) 
at J^ + . (The calculation is modeled after which discussed the relation between the 



two symplectic structures in the absence of internal boundaries within the framework of 
geometrodynamics.) For simplicity of presentation, we will just make assumptions on the 
asymptotic behavior of fields as they are needed in the intermediate stages of the calculation 
and collect our assumptions at the end. 

The canonical symplectic structure ( |6.1U| ) can be obtained by integrating a symplectic 



current to on the partial Cauchy surface M. This 3-form u is given by 

00 = f^Tr[SiA A 5 2 S - 5 2 A A fcE] - ^d[6^ 8 2 ( 2 e) - 5 2 ^ 5 1 ( 2 e)] 

07TG 07TG 

+ —SiA A 5 2 *F - 5 2 A A <5i*F 
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(C.6) 



Note that the exact differential in the expression of uj vanishes at infinity due to the fall-off 
conditions on A. Hence, on integrating over M, it provides just the surface term at the 
horizon in the expression ( p. 10] ) of the symplectic structure. The expression of the 3-form u 
involves 4-dimensional fields. However, when we integrate it over M to obtain the symplectic 
structure, only the pull-backs to M of these fields contribute. 

The main idea behind our calculation can be summarized as follows. When equations of 
motion are satisfied, the 3-form uj is curl-free. Therefore, the integral of du trivially vanishes 
in the 4-dimensional region M. bounded by the isolated horizon A, a partial Cauchy surface 
M and null infinity J? + . Hence, provided all fields remain regular (in a conformal completion 
in which i + is a single point), the integral of uj on M equals the sum of the integrals over 
A and J^ + . Let us first consider the integral over A. Using the isolated horizon boundary 
conditions one can show the sum of the first and last terms can be expressed as an exact 
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differential which is precisely the negative of the one appearing in the second term of 10. 
Thus, the integral of the symplectic current over A vanishes. Hence, the integral of the 
symplectic current over M equals that over J^ + . The former is just the canonical symplectic 
structure (|6.10|) . The idea now is to show that the latter is the radiative symplectic structure 



at J^+. 

Let us therefore evaluate Jy+ to. It is immediate from the falloff condition A ~ 0{\) 
that the exact part (i.e., the second term in the expression) of to does not contribute at J^ + . 
Next, using the conformal invariance of Maxwell's equations, it is fairly straightforward 
to evaluate the electro-magnetic part (i.e. the third term) of this integral. Since A = A 
satisfies Maxwell's equations on (Ai, g a b), the Maxwell potentials A are well-behaved at 
Therefore, we impose the gauge condition (|C.4|) at J? + and evaluate the electro-magnetic 
contribution to the integral of uo over J^ + . It equals precisely the electro- magnetic part of 
the symplectic structure (|C.5|) at J^ + . 

Thus, the non-trivial part of the calculation lies in integrating the first term in the 
symplectic current over ^ + . We will now sketch the main steps. 

The gravitational symplectic structure ( |6.10| ) is expressed in terms of the fields £ and A. 



To compare it with the symplectic structure at , we first need to re-express it in terms 
of the metric g a b and its variations. Since we are assuming the equations of motion, and in 
particular Gauss' law DH = 0, the connection A can be expressed in terms of the soldering 
form a as 

a AB 1 bAA'trr „ B { n 7\ 

K = -7^ V a a 6 A ,. (C.7) 
An arbitrary variation to the soldering form a can be written as 

5a AA ' = \{5g ah a hAA ' + ^ ab a bAA ') (C.8) 

where 5g ab is symmetric and fi ab is antisymmetric. It is easy to check that the above vari- 
ation in a induces a variation 5g ab in the metric. Also, by performing an internal gauge 
transformation on 5a without changing the background field a, fi ab can be set equal to zero. 
This gauge transformation leaves the variation of the metric 5g ab unchanged. Hence, from 
now on, without loss of generality, we assume we are in a gauge in which \i ab has been set 
to zero. This choice of internal gauge will simplify our calculations considerably but is not 
essential since the symplectic structure is gauge invariant. We can now express 5a in terms 
of the conformally rescaled soldering form a a = Qa a as 

^ a AA ' = ^9a b a bAA ' (C.9) 
Now, one would naively expect Q 2 5g a b to be finite at J^ + . However, in the context of 



vacuum general relativity, Geroch and Xanthopoulos |39| have shown that perturbations 
with C°° initial data of compact support have a better behavior: in a suitable gauge, Q 2 5g a b 
in fact vanishes at J? + . Furthermore, 

h ab :=n5g ab (CIO) 
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is well defined and C°° at and satisfies the following conditions: 

QT x n a h ah and Vt- 2 h a h b h ab are and g ab ^nh ab = (C.ll) 

at J? + . Finally, the trace-free part of the field h ab at J? + is precisely the field 5^ ab representing 
the change in the equivalence class of connections at J^ + , i.e. the tangent vector to the phase 
space of radiative modes induced by 5g ab 



5~fab = h ab -^q mn h mn q ab . (C.12) 

We will assume the tangent vectors S under consideration have this asymptotic behavior. 

With this structure at hand, a straightforward but lengthy calculation enables one to 
express the variations of A and £ in terms of fields which are smooth at 

SA C AB = \t^ AB p{V e h cf ) + g cf h ek n k ] 

^ab B = ^^f/ B hb]d 



Using the identity 



^ab^ ^cdAB = 2 (9ac9bd ~ 9ad9bc) ~ ^ abcd (C.14) 



and simple consequences of the Geroch-Xanthopoulos asymptotic behavior ( |C.ll 



one can 



now express the symplectic current in terms of the fields h\ ab and h 2a b- Assuming that at 
least one of the two perturbations, h lab and h 2a b vanishes at i° and i + , we can therefore write 
the gravitational part of the canonical symplectic structure as 

nft.fc) = ^ I + (h lab n c V c h 2cd - h 2ab h c V c h lcd )g ac 9 M 'i, (C.15) 

where the volume form ^ on is given by ( |C.1|) . 

To bring this expression to the same form as appears in ((J. 5) it is necessary to replace 



the fields h ab with their trace-free parts Sjab in the first integral. This will not introduce any 
additional terms because of the properties flC.llp of h ab . Also, since J? + is divergence-free, 
we can replace h c V c h a b with J??nh ab . 

In summary, we have shown that, when the equations of motion hold, both the gravita- 
tional and electro-magnetic parts of the symplectic structure can be rewritten in terms of 
fields living at J^ + . Combining these results, it follows that 

Q{S 1 , 8 2 ) = fi rad (5l ad , ^ ad ), (C.16) 

provided the background and the tangent vectors have certain asymptotic properties. 

To conclude, let us collect the assumptions on the behavior of various fields that were 
necessary to arrive at (|C.16| ). The background solution is assumed to be asymptotically flat at 
spatial and future null infinity and asymptotically Schwarzschild at future time-like infinity. 
In a conformal frame in which J" + is divergence-free, the linearized fields h ab are assumed to 
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satisfy the Geroch-Xanthopoulos conditions ( |C.10|) , (|C.11| ) and the Maxwell potential SA is 
assumed to satisfy (|C.4|) at J? + . Next, at least one of h lab and h 2a b has to vanish at i + and 
at least one of them has to vanish at i°. This last assumption can easily be met in the actual 
application of (|C.16|) in the main text (Section |7.3|) . There, S 2 is the Hamiltonian vector 
field associated with a BMS time translation n a ; S 2 = (^hlab, -^hA). Now, up to numerical 
factors, the total energy radiated across in the (background) space-time is given by the 
integral of squares of these two fields. Hence, it is physically reasonable to restrict oneself 
to space-times in which the two fields go to zero as one approaches i + and i° along J^ + . In 
this case, h 2a b w ih automatically satisfy the last requirement. Finally, for the main result 
( [7.9Q of Section [7.3| to hold, an additional condition must be satisfied: the linearized fields 
h ab and SA satisfying ( |C.10| ), ( |C.11|) and ( |C.4j) should span the tangent space at each point 
of the sector of phase space considered. 

While these assumptions seem plausible, we do not know of general results which will 
ensure that a 'sufficient number' of such background solutions exist or that they will admit 
a 'sufficient number' of linearized fields satisfying our conditions. Indeed, at this stage, 
one does not even have a conclusive proof of existence of a 'sufficient number' of radiating 
solutions which have smooth and complete J? + . 
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